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ABSTRACT 

Difficulties in the subspecific identification of more than 200 specimens of Thomomys um¬ 
brinus from the State of Tlaxcala led us to examine the systematic status of the eight subspecies 
recognized within the Mexican Transvolcanic Belt. We gathered representative samples of 
those subspecies (> 50) during three years of fieldwork. These were analyzed together with 
samples from other populations not previously assigned to subspecies, as well as voucher 
specimens from several collections. Thirty-seven variables, obtained from 1145 specimens, 
were analyzed morphometrically using univariate and multivariate statistics. Within-population 
variation was highly influenced by sexual dimorphism, and especially by ontogenetic variation; 
therefore, these two factors should be considered carefully when making systematic and tax¬ 
onomic decisions. Our analysis of geographic variation determined only two recognizable 
subspecies of Thomomys umbrinus within the Mexican Transvolcanic Belt: Thomomys umbri¬ 
nus umbrinus (Richardson, 1829) and Thomomys umbrinus pullus Hall and Villa, 1948. 


RESUMEN 

Dificultades en la identificacion subespecffica de mas de 200 especfmenes de Thomomys 
umbrinus del estado de Tlaxcala, nos llevaron a examinar el estado sistematico de las ocho 
subespecies reconocidas en el Eje Volcanico Transversal. Obtuvimos muestras representativas 
de esas subespecies (> 50) durante tres anos de trabajo de campo. Estas fueron analizadas 
junto con muestras de otras poblaciones no asignadas previamente a las subespecies y ejern- 

1 Departamento de Biologfa, Universidad Autonoma Metropolitana Unidad Iztapalapa. Apdo. Postal 55-535. 09340 
Mexico D. F. 

2 Research Associate, Division of Vertebrate Zoology (Mammalogy), American Museum of Natural History. 


Copyright © American Museum of Natural History 2000 


ISSN 0003-0082 / Price $4.70 



2 


AMERICAN MUSEUM NOVITATES 


NO. 3297 


plares de referenda de varias colecdones. Treinta y siete variables fueron obtenidas de 1145 
ejemplares y fueron analizadas morfometricamente utilizando estadlstica uni y multivariada. 
La variacion intrapoblacional resulto altamente infiuenciada por el dimorfismo sexual, pero 
especialmente por la variacion ontogenetica; por ende, estos dos factores deben ser conside- 
rados cuidadosamente cuando se tomen decisiones sistematicas y taxonomicas. Nuestro analisis 
de la variacion geografica determino que solo pueden reconocerse dos subespecies de Thom¬ 
omys umbrinus en el Eje Volcanico Transversal: Thomomys umbrinus umbrinus (Richardson, 
1829) y Thomomys umbrinus pullus Hall y Villa, 1948. 


INTRODUCTION 

A large and highly heterogeneous sample 
of Thomomys umbrinus from Tlaxcala in¬ 
cluded the whole variety of diagnostic char¬ 
acters for the eight subspecies acknowledged 
in the MTB, Mexican Transvolcanic Belt, 
(Nelson and Goldman, 1934; Hall and Villa, 
1948). In fact, the Tlaxcalan specimens could 
be assigned to any of the following subspe¬ 
cies: Thomomys umbrinus albigularis from 
Hidalgo; Thomomys umbrinus martinensis 
and Thomomys umbrinus orizabae from 
Puebla; or Thomomys umbrinus vulcanius 
from Mexico (Hall, 1981). 

Examination of most of the holotypes and 
topotypes of the subspecies from the MTB 
housed at National Museum of Natural His¬ 
tory (USNM) and the American Museum of 
Natural History (AMNH), made clear the 
qualitativeness of the diagnostic characters. 
For example, Thomomys umbrinus, vulcan¬ 
ius, T. u. tolucae, T. u. martinensis, T. u. um¬ 
brinus, T. u. albigularis, and T. a. peregrinus 
are similar both in shade and in cranial char¬ 
acteristics. In fact, some T. u. albigularis 
lack the allegedly diagnostic white spot on 
their throat (Nelson and Goldman, 1934), but 
some T. u. peregrinus and all of those of T. 
u. vulcanius do have it. In T. u. orizabae, the 
holotype and most of the topotypes differ 
from the preceding by their melanic shade; 
however, melanic individuals were also 
found among Tlaxcalan specimens as well as 
in T. u. pullus collected in Michoacan. 

Those qualitative characters, together with 
the apparently isolated distribution of the 
populations along the MTB, were previously 
used in designating seven of the eight sub¬ 
species in the MTB (Merriam, 1893; Bailey, 
1915; Nelson and Goldman, 1934). Since 
only small sample sizes (n < 15) were ex¬ 
amined, investigators emphasized a geo¬ 


graphic but not necessarily a genetic isola¬ 
tion nor a wider morphological variation. 

Therefore, we decided to allocate the Tlax¬ 
calan specimens according to the morpho¬ 
metric variation of the subspecies within the 
MTB, expecting that our results would mo¬ 
tivate a thoughtful study of variation in this 
taxon throughout its range. Such scope 
would enable greater comprehension of the 
systematics and taxonomy of Thomomys um¬ 
brinus in the MTB. It might also contribute 
to a better understanding of the role of sub¬ 
species as a means of recognizing geographic 
variants of the species (Mayr, 1969). This is 
possible if the subspecies is regarded as con- 
specific populations sharing their genetic 
pool through gene flow (Smith and Patton, 
1988; Patton and Smith, 1989), on one hand, 
but which have a particular morphometric 
variation because they are adjusted to the 
habitat conditions in their geographic area 
(Grinnell, 1935; Yates and Schmidly, 1977; 
Honeycutt and Schmidly, 1979; Patton and 
Smith, 1990). 

After three years of fieldwork, we gathered 
large sample sizes (n > 50) along the MTB. 
In addition, voucher specimens from the 
MTB were examined in several collections. 
Qualitative examination of collected and 
voucher specimens showed that persistence 
of morphological features was not consistent 
within a particular locality; on the contrary, 
the variation was generally noticeable. The 
only exception was the distinctly small-sized 
specimens of T. u. pullus. Moreover, in ad¬ 
dition to the overlapping morphological var¬ 
iation among the samples from MTB, stan¬ 
dard karyotypes and allozymic data showed 
no distinctions among some of them (Hafner 
et al., 1987). That is, T. u. pullus (Patzcuaro, 
Michoacan), T. u. tolucae (Toluca, Mexico), 
T. u. vulcanius (Amecameca, Estado de Mex¬ 
ico), and T. u. umbrinus (Boca de Monte, 
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Veracruz) belong to the same chromosomic 
(NF = 78) and allozymic groups (Hafner et 
al., 1987). These data together certainly do 
not suggest distinct species, and perhaps not 
even separate subspecies. 

Therefore, the objective here was to ex¬ 
amine geographic variation among the pop¬ 
ulations of Thomomys umbrinus in the Mex¬ 
ican Transvolcanic Belt. Morphometric anal¬ 
yses include the eight previously recognized 
subspecies and populations not previously al¬ 
located from the states of Mexico, Michoa- 
can, and Tlaxcala. 

Diagnostic Characters, Taxonomy, and 
Distribution 
of Thomomys umbrinus 

Thomomys umbrinus is a rodent species of 
fossorial habits belonging to the Geomyinae 
within the Geomyidae. The genus is charac¬ 
terized (Hall, 1981) by its small size and by 
the rather delicate structure of its hands with 
fine claws, as opposed to the robust hands 
with strong claws of other genera ( Geomys 
and Cratogeomys). The genus Thomomys 
does not show the basitemporal fossa located 
between the lower third molar (m3) and the 
lingual face of the coronoid process is rather 
a slight depression in some specimens in the 
United States. Thomomys has the front sur¬ 
face of the upper incisors smooth, but with a 
very fine groove toward the inner limit, 
while the other genera of the family have ei¬ 
ther one ( Cratogeomys , Pappogeomys and 
Zygogeomys) or two well-defined ( Orthogeo- 
mys ) grooves. In addition, the third upper 
molar (M3) is elliptic, monoprismatic, and 
similar to the other two upper molars (Ml, 
M2) in Thomomys. Conversely, in other gen¬ 
era the M3 is rather biprismatic and differs 
from the Ml and M2 (Hah, 1981). 

Within the genus, systematic and taxo¬ 
nomic status of Thomomys umbrinus has un¬ 
dergone numerous changes (Hah and Kelson, 
1959; Anderson, 1966; Hoffmeister, 1969; 
Patton and Dingman, 1968; Patton, 1973; 
Hall, 1981; Ramirez-Pulido et al., 1996). 
Such changes include: (1) synonymy with 
other congeneric species (Bailey, 1915; Hah 
and Kelson, 1959; Hall, 1981) and (2) rec¬ 
ognition of Thomomys bottae and T. town- 
sendii as separate species (Anderson, 1966; 


Hoffmeister, 1969; Patton and Dingman, 
1968; Patton, 1973; Patton and Smith, 1990). 
Furthermore, changes include the subsequent 
reorganization of the involved subspecies. 

Before conclusive evidence was gathered 
(Patton and Dingman, 1968; Patton, 1973; 
Patton and Smith, 1990) to formally consider 
them as separate species (Patton, 1993), Hah 
(1981) combined Thomomys umbrinus with 
Thomomys bottae and included 229 subspe¬ 
cies. Therefore, the distribution of Thomomys 
umbrinus started north of the United States 
and extended along the Californian mountain 
ridges and the Sierra Madre Occidental and 
Sierra Madre Oriental, all the way down to 
the MTB. 

The distribution is discontinuous, particu¬ 
larly in the states of the Mexican Plateau and 
the MTB. The presence of arid zones and, 
probably, the presence of some species of 
Cratogeomys (Hah, 1981) have contributed 
to this pattern. Indeed, Nelson and Goldman 
(1934) noted that Cratogeomys merriami 
merriami restricted the distribution of Thom¬ 
omys at 11,500 ft in the Popocatepetl Vol¬ 
cano, and that Cratogeomys tylorhinus plan- 
iceps did the same at 8600 ft in the Nevado 
de Toluca Volcano. Observations by Nelson 
and Goldman (1934) were confirmed during 
fieldwork, especially in the Valley of Mexi¬ 
co, where Cratogeomys occurs in many lo¬ 
calities in which Thomomys umbrinus is no 
longer found. 

Simpson (1961) strongly criticized the 
large number of subspecies allocated to 
Thomomys umbrinus (Smith and Patton, 
1988), and in so doing, he placed this taxon 
in the midst of a debate on the validity of 
the subspecies (Grinnell, 1935; Davis, 1938; 
Wilson and Brown, 1953; Brown and Wil¬ 
son, 1954; Simpson, 1961; Choate and Wil¬ 
liams, 1978; Smith and Patton, 1988; Patton 
and Smith, 1990). In fact, studies by Patton 
et al. (Patton, 1973; Patton and Feder, 1978; 
Hafner et al., 1987) have raised new ques¬ 
tions about the taxonomic and systematic 
population status of this species. Hafner et 
al. (1987) found three chromosomic groups 
within Thomomys umbrinus that show genic 
differentiation at the allozymic level, sug¬ 
gesting more than one species of pocket go¬ 
pher as currently recognized. 

After the separation of Thomomys umbri- 
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nus and Thomomys bottae, 47 subspecies 
from the 78 occurring in Mexico (Hall, 1981) 
have been allocated to the latter (Patton and 
Smith, 1990; Patton, 1993; Rarmrez-Pulido 
et al., 1996). Therefore, T. umbrinus is rep¬ 
resented by the remaining 31 subspecies 
which occur along the Sierra Madre Occi¬ 
dental from southern Sonora downward to 
Durango and Nayarit. In the Mexican Pla¬ 
teau, the species is distributed in patchy ar¬ 
rangements in Chihuahua and Coahuila in 
the north but in more isolated populations in 
Zacatecas, Aguascalientes, San Luis Potosf, 
and Guanajuato in the center. In the Sierra 
Madre Oriental, T. umbrinus extends from 
mid-Coahuila to Nuevo Leon. The species 
was reported to occur in the MTB within iso¬ 
lated populations in Michoacan, Mexico, 
Distrito Federal, Morelos, Tlaxcala, Puebla, 
Hidalgo, and Veracruz (Hall, 1981). 

The eight subspecies occurring in the 
MTB (Hall, 1981) are Thomomys umbrinus 
umbrinus (Richardson) in Veracruz; Thomo¬ 
mys umbrinus orizabae Merriam in Puebla; 
Thomomys umbrinus peregrinus Merriam in 
Distrito Federal, Hidalgo, Mexico, and Mo¬ 
relos; Thomomys umbrinus albigularis Nel¬ 
son and Goldman in Hidalgo, Puebla, and 
Veracruz; Thomomys umbrinus martinensis 
in Puebla; Thomomys umbrinus tolucae Nel¬ 
son and Goldman in the state of Mexico; 
Thomomys umbrinus vulcanius Nelson and 
Goldman in Mexico, Morelos, and Puebla; 
and Thomomys umbrinus pullus Hall and 
Villa in Michoacan. 

Study Area 

The Mexican Transvolcanic Belt (MTB) 
includes the highest regions of Mexico (2000 
to 5650 m), located between 19° and 20°N 
and 96° and 105°W, and crossing the country 
from ocean to ocean. From west to east, the 
volcanoes related to this study are Nevado 
de Toluca (4624 m), Popocatepetl (5452 m), 
Iztaccihuatl (5460 m). La Malinche (4461 
nr), and Pico de Orizaba (5700 m). The re¬ 
gion is 130 km deep and more than 1000 km 
long. It is bounded by the Sierra Madre Ori¬ 
ental in the east, by the Sierra Madre del Sur 
in the south, and by the Sierra Madre Occi¬ 
dental in the west. There are about 1000 vol¬ 
canoes from the Quaternary Age of less than 


1 million years (Fa and Morales, 1991, De- 
mant, 1978; Nixon et al., 1987; Pasquare et 
al., 1987a, 1987b) in the MTB. 

Specimens were collected in the second 
through fifth districts recognized by Demant 
(1978) in the MTB. The second district in 
the east shows numerous monogenetic cones 
and, except the Tancftaro Hill in Michoacan, 
there are no stratified volcanoes. Four large 
stratified volcanoes are located from the third 
to the fifth districts at the center of the MTB 
(Nevado de Toluca, Popocatepetl, Iztaccf- 
huatl, and La Malinche) and oriented in a 
north-south direction within the lake depres¬ 
sion corresponding to the Sierra Nevada. 
Other small sierras, among which Chichi- 
nautzin is found, run in a northeast-southwest 
direction, on fissures closer to the southern¬ 
most portion of the Valley of Mexico. In the 
fifth district, the Cofre de Perote-Pico de Ori¬ 
zaba volcano chain, which is a recent lava 
cone on an older volcano, runs eastward. 

The sampled zones have a temperate cli¬ 
mate with rainfall at the end of spring and 
fall, and the northern portions are less humid 
than their southern counterparts. The annual 
average temperature is generally about 18°C, 
but it decreases as altitude increases (Garcia, 
1981). The predominant vegetation in the 
studied areas (Rzedowski, 1978) is oak forest 
(Quercus sp.), pine-oak ( Quercus-Pinus sp.) 
and pine-spruce ( Pinus-Abies religiosa) as 
well as pine forest with fodder grass ( Muhl- 
enbergia sp. and Festuca sp.). Because the 
parent rock is volcanic, the ground is sandy 
near the volcanic uplifts and more clayey in 
the lower regions. 

METHODS 

Collecting and Examining Specimens 

Previously, specimens from the study area 
were less than 15 in most collections and no 
extensive morphometric work had been con¬ 
ducted. Therefore, we decided to collect in¬ 
tensively and to examine voucher specimens 
in several mammal collections. After three 
years of exhaustive fieldwork, more than 50 
representative samples of most of the rec¬ 
ognized subspecies were obtained from the 
type locality and neighboring areas. The two 
exceptions were T. u. martinensis (n = 47) 
and T. u. umbrinus (n = 5). Specimens were 
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also collected from unsampled areas between 
populations of the recognized subspecies in 
the states of of Mexico, Tlaxcala, and Mi- 
choacan. Such samples suggest gene flow 
among all populations of Thomomys umbri¬ 
nus in the MTB. In total, 857 specimens 
were gathered from 67 localities in the states 
of Hidalgo, Mexico, Michoacan, Morelos, 
Puebla, Tlaxcala, and Veracruz. The collect¬ 
ed specimens were prepared as skins and 
skulls (Hall, 1981; Rarmrez-Pulido et al., 
1989) and they are housed at Coleccion de 
Mamiferos de la Universidad Autonoma Me- 
tropolitana, Unidad Iztapalapa (UAMI). The 
skulls were cleaned with dermestid beetles 
(Hall and Russell, 1933; Rarmrez-Pulido et 
al., 1989). 

In addition, 288 specimens from nine lo¬ 
calities, were examined in other collections: 
Coleccion de Mamiferos del Laboratorio de 
Cordados Terrestres, Escuela Nacional de 
Ciencias Biologicas, Instituto Politecnico 
Nacional (CB) at Mexico City; American 
Museum of Natural History (AMNH) at New 
York; National Museum of Natural History, 
Smithsonian Institution (USNM) at Washing¬ 
ton, D.C.; Texas Cooperative Wildlife Col¬ 
lections (TCWC) of the Texas A&M Uni¬ 
versity in College Station; Museum of Nat¬ 
ural History (KU) of the University of Kan¬ 
sas, Lawrence. Among specimens in other 
museums, most holotypes (USNM) and sev¬ 
eral topotypes (AMNH, CB, KU, UAMI) of 
the mentioned species were examined. Ex¬ 
ceptions include the holotype of Thomomys 
umbrinus pullus, at the University of Mich¬ 
igan, and that of Thomomys umbrinus um¬ 
brinus, at the Natural History Museum, Lon¬ 
don. 

We examined 1145 specimens of Thomo¬ 
mys umbrinus from 76 localities in the six 
mentioned states. However, the number of 
specimens included varied with the type of 
statistical analysis. In addition, some very 
young specimens or damaged skulls were 
eliminated. The numbers of specimens in 
each calculation is mentioned in the methods 
or in the tables for each study. 

Sex of most specimens was gathered from 
museum catalogs and labels. In the few cases 
in which sex was unrecorded, the skins were 
examined for evidence such as well-devel¬ 
oped scrotum or nipples. The pelvic structure 


was also checked for sexual differences (i.e., 
a larger pubic foramen and a wider simphysis 
in females); however, if the postcranial skel¬ 
eton was not available, the specimen was ex¬ 
cluded. 

Hair color of holotypes, topotypes, and 
other specimens (n = 76, 55 2; 21c?) was 
reviewed. We looked for qualitative taxo¬ 
nomic characteristics that might distinguish 
the subspecies. Likewise, in order to obtain 
data about the habitat and the exact capture 
localities, the original collectors’ notes on the 
specimens were reviewed. Among others, 
notes of W. Dalquest in 1946 (KU), S. An¬ 
derson in 1952 (KU), and E. W. Nelson in 
1892 (USNM), were examined. 

Measurements 

Lour conventional external measurements 
(Hall, 1981) together with 33 skull measure¬ 
ments (fig. 1) were recorded from each spec¬ 
imen. The external measurements are the to¬ 
tal length of the body (TOL); length of the 
vertebral tail (LVT); length of the right hind 
foot (HLT); and length of the ear (LOE). 
Measurements were taken from the skin la¬ 
bels of voucher specimens. 

Some skull measurements were taken from 
the literature (Thaeler, 1968; Heaney and 
Timm, 1983) and others were suggested by 
Dr. Michael Smolen (Texas A&M Universi¬ 
ty, Nov. 1989). However, most measure¬ 
ments were selected during direct examina¬ 
tion of the specimens by one of us (ACC). 
Cranial measurements (fig. 1) were taken 
with calipers and recorded to the nearest 0.01 
mm. The following list contains the names 
and abbreviations of measurements (numbers 
in parentheses correspond to fig. 1). 

Braincase and Rostrum 

GLS (1-1') greatest length of skull 

LNA (1-2) length of nasals 

LMX (2'—2") length of maxilla 

LFO (3-3') length of frontal 

LZA (4—4') length of zygomatic arch 

LBC (4—5) length of braincase 

ICO (6-6') interorbital constriction 

WRO (7-7') width of rostrum 

AZB (8-8') anterior zygomatic breadth 

MZB (9—9') medial zygomatic breadth 

PZB (10—10') posterior zygomatic breadth 



6 


AMERICAN MUSEUM NOVITATES 


NO. 3297 


SBP (11—11') skull breadth across parietals 
(i.e., braincase) 

LDI (12—12') length of diastema 

LMT (12-13) length of molar toothrow 
LPF (14—14') length from the pterygoid fossa 
to the basioccipital 

BAM (15-15') breadth across molars (Ml-Ml) 
BAB (16-16') breadth across bullae 

BOP (17—17') breadth across basioccipital pro¬ 

cesses 

HRO (18-18') height of rostrum 

HMC (19-19') height of maxillary crests 

MHS (20-20') maximum height of skull 

HBP (21—21') height from the bulla to the pa¬ 
rietal 

Mandible 

BCP (22—22') breadth across coronoid pro¬ 
cesses 

BJT (23—23') breadth across jaw teeth (m2- 
m2) 

AWD (24-24') anterior width of dentary 
GWM (25-25') greatest width of mandible 
MWD (26-26') medial width of dentary 
GLD (27-27') greatest length of dentary 
LCP (28—28') length to coronoid process 
BML (28—29) basal mandible length 
DHI (27—30) dentary height from basis of in¬ 
cisive alveolum 

AHM (28-31) anterior height of mandible 
PHM (29-28') posterior height of mandible 

Sampled Localities 

Sampling was carried out within: (a) lo¬ 
calities represented by few specimens (n < 
50); (b) localities in potential corridors for 
gene flow among previously recognized sub- 
species; (c) localities still occupied by Thom¬ 
omys umbrinus', or (d) topotypic localities of 
subspecies not represented in Mexican col¬ 
lections (i.e., T. n. albigulciris, T. u. marti- 
nensis, T. u. orizabae, T. u. pulliis, T. u. to- 
lucae, and T. u. umbrinus). 

In the field, we noticed that T. u. albigu- 
laris, T. u. pulliis, and T. u. martinensis were 
becoming scarce at the type localities. This 
rarity results from habitat modifications and 
from displacement by Cratogeomys merria- 
mi. In fact, only one topotypical female of T. 
u. albigularis from the National Park “El 
Chico,” at the Sierra de Pachuca, Hidalgo, 
was collected and most of the 58 inividuals 
of T. u. pullus were collected at Los Tanques, 
Michoacan. We collected no specimens of T. 
u. martinensis at San Martin Texmelucan, 


State of Mexico, during the three years of 
fieldwork. Fortunately, 40 individuals refer¬ 
able to T. u. martinensis were located at the 
USNM. 

Nelson and Goldman (1934) described 
Thomomys umbrinus orizabae and, especial¬ 
ly, T. u. tolucae and T. u. vulcanius with im¬ 
precise localities. Therefore, our sampling 
was conducted as close as possible to the lo¬ 
calities described in the fieldnotes. Accord¬ 
ingly, specimens of T. u. orizabae were ex¬ 
amined from Miguel Hidalgo y Costilla, 
Puebla. In the state of Mexico, we collected 
specimens of T. u. tolucae at several locali¬ 
ties near the Nevado de Toluca Volcano, and 
the same was done with T. u. vulcanius at 
the Paso de Cortes. 

Sampled localities (SL) are mapped in fig¬ 
ure 2 and are detailed under Examined Spec¬ 
imens. Samples of Thomomys umbrinus col¬ 
lected in zones not previously recorded, but 
which were in the gaps between two recog¬ 
nized subspecies, are considered to be vouch¬ 
ers for corridors of genetic continuity be¬ 
tween the eight subspecies inhabiting the 
MTB (fig. 2). 

Nongeographic Variation 

Individual variation, ontogenetic variation, 
and sexual dimorphism were examined in the 
population of the Sierra de Tlaxco, Tlaxcala 
(n = 356). We considered four age groups (1 
= young; 2 = subadults; 3 = adults; 4 = old 
specimens) among 214 females and 142 
males according to the expansion of the zy¬ 
gomatic arches and the elongation of the ros¬ 
trum relative to the braincase (Castro-Cam- 
pillo et al., 1993). 

Individual variation in each external and 
cranial measurement was analyzed by obtain¬ 
ing the standard statistics (mean, x; standard 
error, SE; minimum value, MIN; maximum 
value MAX; coefficient of variation, CV) 
with the routine PROC UNIVARIATE SAS 
(SAS, 1985). These analyses also served to 
corroborate the normal distribution of the 
data and helped in finding outlier values. 

We determined ontogenetic or age varia¬ 
tion and the presence of secondary sexual di¬ 
morphism for each taxonomic character, 
through a one-way analysis of variance with 
ANOVA routine. When significant differen¬ 
ces (p < 0.05) between ages or sexes were 
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Fig. 1. Cranial measurements (33) taken in 1145 Thomomys umbrinus from 72 localities in the 
Mexican Transvolcanic Belt. A. a: dorsal view; B, b: ventral view; C, c: lateral view. See pp. 13 - 15 
for explanations of measurements. 
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TABLE 1 

The 19 Group Localities Sampled in the Mexican Transvolcanic Belt 

Abbreviations (Code) refer to the 17 taxa of Thomomys umbrinus included in this study, 
n = sample size. D.E = Distrito Federal. See map in figure 2. 


Code 

Taxon 

Group locality 

State 

n 

P 

T. u. pullus 

Patzcuaro 

Michoacan 

58 

N 

T. umbrinus ssp. 

Cerro San Andres 

Michoacan 

27 

L 

T. umbrinus ssp. 

Lengua de Vaca 

Mexico 

25 

B 

T. umbrinus ssp. 

Amanalco de Becerra 

Michoacan, Mexico 

67 

C 

T. umbrinus ssp. 

Sierra de las Cruces 

Mexico 

42 

O 

T. u. tolucae 

Volcan Nevado de Toluca 

Mexico 

102 

E 

T. u. peregrinus 

Salazar 

Mexico, D.E, Morelos 

72 

V 

T. u. vulcanius 

Volcan Popocatepetl 

Mexico, Puebla 

88 

M 

T. u. martinensis 

San Martin Texmelucan 

Puebla 

47 

A 

T. u. albigularis 

Sierra de Pachuca, El Chico 

Hidalgo, Puebla 

82 

X 

T. umbrinus ssp. 

Sierra de Tlaxco 

Tlaxcala, Puebla 

407 

I 

T. umbrinus ssp. 

Volcan La Malinche 

Tlaxcala, Puebla 

92 

R 

T. u. orizabae 

Volcan Pico de Orizaba 

Puebla 

62 

S 

T. u. umbrinus 

Boca de Monte 

Puebla, Veracruz 

5 

T 

T. umbrinus ssp. 

Totoltepec 

Puebla 

5 

D 

T. umbrinus ssp. 

La Soledad 

Tlaxcala 

12 

Y 

T. umbrinus ssp. 

Oyameles 

Puebla 

8 


found, Duncan’s Test for Multiple Means or 
Tukey’s Test for the Student-Type Range 
(PROC GLM) were used, respectively (Ott, 
1984), in order to determine subgroups 
among the means. 

We considered age variation separately for 
males and females and, at the beginning, all 
of the specimens were included. However, 
since young (age class 1) became separated 
from the others in 86.65% of all features, we 
eliminated class 1. Both analyses are pre¬ 
sented, but emphasis is given from subadults 
on (age classes 2—4). Accordingly, for sec¬ 
ondary sexual dimorphism, analyses includ¬ 
ed each age group separately, and the young 
specimens were excluded. 

We ran two-way analysis of variance fol¬ 
lowed by Duncan’s test (PROC GLM) to de¬ 
termine the effect of sex, age, and the inter¬ 
action between sex and age on the variance 
in the examined sample (Hollander, 1990; 
Stangl et al., 1991). In Duncan’s tests on sex 
and age, only the specimens aged 2 and up 
were considered. Finally, in order to under¬ 
stand the role played by sex, age, and their 
interaction within total nongeographic vari¬ 
ation in Thomomys umbrinus, we analyzed 
the components of the variance (Straney, 
1978; Leamy, 1983; Patton and Smith, 1990) 
using the VARCOMP routine. 


Geographic Variation 

We examined the morphometric similari¬ 
ties among the populations of Thomomys um¬ 
brinus from the MTB by arranging the sam¬ 
ple localities (n = 76) into 17 group localities 
(GL, table 1, fig. 2). Both nearness and eco¬ 
logical homogeneity were considered in at¬ 
tempting to select groups with no geographic 
or ecological barriers that would restrict gene 
flow within each group (Patton and Smith, 
1990). 

Based on the analysis of intrapopulation 
variation, we calculated the nature and mag¬ 
nitude of geographic variation of T. umbri¬ 
nus in the MTB separately for both sexes and 
age classes (2—4). 

Statistical programs exclude individuals 
lacking information in any of the variables; 
therefore, to include the greatest number of 
specimens and GLs, a regression analysis 
(PROC GLM, SAS) was developed to esti¬ 
mate the missing values in individuals. Spec¬ 
imens lacking more than two measurements 
were excluded from the analysis (n = 30). 
Likewise, some measurements were elimi¬ 
nated because they were missing in several 
specimens (BCP, BJT, AWD, GWM, MWD); 
were too variable (TOL, LVT, HFT, LOE; 
BOP; and BML); contributed no information 
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to differentiate among the GLs; or were high¬ 
ly correlated to other variables like MWZ 
and PWZ (R = 0.8, Castro-Campillo et al., 
1993). 

Geographic variation was examined by 
univariate and multivariate statistical analy¬ 
ses. Standard statistics were calculated for 
each variable with UNIVARIATE (SAS) and 
then a one-way ANOVA (GLM, SAS) was 
run to disclose significant (p < 0.05) inter¬ 
population differences among means of the 
GLs. This analysis was completed by Dun¬ 
can’s Test for Multiple Means (SAS) to find 
population subgroups with no significant dif¬ 
ferences. The ANOVA was made on the GLs 
having four or more specimens. 

Principal Components Analysis (PCA) 
with the PRINCOMP (SAS) routine on the 
measurements of each specimen determined 
the trends in the variation pattern. Results are 
summarized in a bidimensional graphic of 
the first (PC I) and third (PC III) principal 
components where the distribution of each 
GL is represented by polygons with their re¬ 
spective centroids. We present PC III instead 
of PC II in this plot for the sake of clarity, 
since in the plot of PC I vs PC II, there was 
an overall massive overlapping of the poly¬ 
gons. We developed another PCA on the 
population measurements of the GLs and for 
the tridimensional representation of this anal¬ 
ysis, we used the statistical package Numer¬ 
ical Taxonomy System NT-SYS (Rohlf and 
Kishpaugh, 1972; Rohlf, 1988). A minimum 
spanning distance tree (MST), calculated 
from the correlation matrices, was projected 
onto the first three PCs (Gower and Ross, 
1969; Rohlf, 1970, 1973, 1975). In that anal¬ 
ysis, as in the rest of the multivariate anal¬ 
yses, all specimens were included. 

Cluster analyses were run on the sample 
means of all the variables in each GL (Sneath 
and Sokal, 1973), using NT-SYS program. 
Group populations with a greater phenetic 
similarity were obtained, using the average 
taxonomic distance as a similarity measure¬ 
ment and the unweighted pair-group method 
based on arithmetic averages (UPGMA). We 
followed this procedure for all age groups 
with sexes separate, as well as in 12 localities 
where subadult (2) and adult (3) specimens 
were collected. The cophenetic coefficient (r) 
of each phenogram was calculated (Lapointe 


and Legendre, 1992) and to compare the re¬ 
sults, Mantel’s tests were developed on the 
correlation and similarity matrices (Mantel, 
1967; Rohlf and Fisher, 1968; Sokal, 1979; 
Smouse et al., 1986). 

Furthermore, a classification discriminant 
analysis (DISCRIM, SAS) was used to de¬ 
termine the number of specimens correctly 
assigned according to the acknowledged des¬ 
ignation or to the GL of origin (Castro-Cam¬ 
pillo et al., 1993). This analysis also deter¬ 
mined the degree of overlap among the GLs 
because of misclassified specimens (Castro- 
Campillo et al., 1993). 

Finally, a brief analysis related pelage col¬ 
oration to original descriptions (Nelson and 
Goldman, 1934; Hall and Villa, 1948). For 
this analysis, only the type specimens and the 
material deposited at UAMI were reviewed. 
Pelage shade is described according to 
Smithe (1974a, 1974b) with the names and 
color numbers in parentheses after the de¬ 
scription. 

Specimens Examined 

Group localities of previously recognized 
subspecies and other populations of Thomo¬ 
mys umbrinus examined along the MTB are 
listed from west to east according to table 1 
and figure 2. Taxonomic names are followed 
by the abbreviated (table 1) and the spelled- 
out name corresponding to the group locali¬ 
ties (GLs) mentioned in tables, figures, and 
text. Next, the total number of specimens ex¬ 
amined is indicated in parentheses. Within 
the sampled localities (SL), the type locality 
is indicated by an asterisk and if the holotype 
was examined, its sex, age, catalog number, 
and the abbreviation for its institutional lo¬ 
cation are indicated also after the name of 
the GL. The SLs are preceded by their as¬ 
signed number according to figure 2 and fol¬ 
lowed by the number of examined males (c3 ) 
and females ( 2 ), and housing museum in pa¬ 
rentheses. 

Thomomys umbrinus pullus. P = Pdtzcuaro 
(45). Michoacan: 2.-7 km SE Patzcuaro, 
2380 m * (8 2 5(3 CB); 3. - Los Tanques, 
2450 m (212 11 <3 UAMI). 

Thomomys umbrinus ssp. N = Cerro de San 
Andres (26). Michoacan: 4. — Laguna 
Verde, 12.5 km N, 18 km W Cd. Hidalgo, 
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2790 m (1 5 1 c? UAMI); 5. - Los Azufres, 
2790 m (15 UAMI); 5. - Los Azufres, 
2850 m (6 2 5 c? UAMI); 6. - Pucuato, 8 
km S, 15 km W, Ciudad Hidalgo, 2400 m 
(5 2 7 c? UAMI). 

Thomomys umbrinus ssp. L = Lengua de 
Vaca (25). Michoacan: 7. — Cerro Prieto, 
8 km S, 3 km W Tlalpujahua, 2780 m (1 2 
UAMI). Mexico: 8. — Laguna del Carmen, 
2800 m (7 2 2c? UAMI); 9. - Cuesta del 
Carmen, 3 km N, 1 km W, Lengua de 
Vaca, 2260 m (lc? CB as peregrinus); 10. 

- 7.5 km S, 2 km E Palizada, 2000 m (1 2 
1 <? UAMI); 11. - 8 km S Lengua de Vaca, 
2480 m (6 2 5 c? UAMI); 12. - 4.5 km S, 
10 km W Villa Victoria, 2560 m (1 2 CB). 

Thomomys umbrinus ssp. B = Amanalco de 
Becerra (59). Michoacan: 13. — 8 km E 
Zitacuaro, 2210 m (25 2 14c? UAMI); 13. 

— 9 km E Zitacuaro, 2300 m (2 2 1 S 
UAMI). Mexico: 14. - 12 km NW Ama¬ 
nalco de Becerra, 2720 m (9 2 lc? UAMI); 
15. - 10 mi N, 6 mi E Valle de Bravo, 
7640 ft (22 3c? KU); 16. - 1.5 mi S Valle 
de Bravo 6050 ft (7 2 5 c? KU). 

Thomomys umbrinus ssp. C = Sierra de las 
Cruces (42). Mexico: 17. — 2.7 km N, 9 
km W Villa del Carbon, 2670 m (1 2 CB 
as peregrinus)-, 18. — Monte de Pena, 2630 
m (212 11c? UAMI); 19. - 9 km N Jilot- 
zingo, 3260 m (3 2 UAMI); 20. — Jilotzin- 
go, 3260 m (3 2 3 c? UAMI). 

Thomomys umbrinus tolucae. O = Nevado 
de Toluca (101). Holotype: Male, age 4 
(55908 USNM).- Mexico: 21. - 3 km SE 
Zinacantepec, 3000 m (28 2 15c? UAMI); 
22. — 5 km S, San Juan de las Huertas, 
3150 m (2 2 7 c? CB as peregrinus)-, 23. — 
Nevado de Toluca, 16 mi SSW Toluca (7 2 
2c? USNM); 23. - Nevado de Toluca, 16 
mi SSW Toluca (lc? CB); 24. - Nevado de 
Toluca, 4 mi S Ralces Toluca (3 2 4 c? 
USNM); 25. - Nevado de Toluca (NW 
slope), 11200 ft (5 2 2c? KU); 25. - N 
slope Nevado de Toluca, 9500 ft * (12 2 
3 c? USNM); 26. - La Saavana, 4 km N 
Ralces, 2500 m (4 2 lc? UAMI); 27. - 
Nevado de Toluca, 6 km S, 5 km W Ra¬ 
lces, 4080 m (4 2 UAMI); 28. - Nevado 
de Toluca, 13 mi S Toluca (lc? USNM). 

Thomomys umbrinus peregrinus. E = Sala¬ 
zar (69). Holotype: Female, age 3 (50130 
USNM).- Mexico: 29. - 7.5 km E Tenango 


de Arista, 2600 m (5 2 1c? CB); 30. - 1 
mi WSW Salazar, 9500 ft (12 lc? KU); 
30. - 1 mi W Salazar, 9850 ft (3 2 4 c? KU); 

30. - 1 mi W. Salazar, 9800 ft (2 2 KU); 

31. - Salazar, 10300 ft * (122 2c? USNM); 
31. - Salazar, 9000 ft (12 USNM); 31. - 
Valle de Toluca, 9000 ft (12 USNM); 31. 

- Salazar, 3050 m (4 2 2c? CB); 32. - La 
Marquesa, 2950 m (9 2 5c? CB). Morelos: 
33. - 5 km N Tres Cumbres, 10,200 ft (4 2 
8<? TCWC); 35. - Cicitepec, 15 km S, 5 
km E Milpa Alta (12 CB). Distrito Fe¬ 
deral: 34. - 2 mi SSW Parres, 9800 ft (2 2 
1 <? KU as peregrinus). 

Thomomys umbrinus vulcanius. V = Popo¬ 
catepetl Volcanoe (88). Holotype: Female, 
age 2 (51885 USNM).- Mexico: 36. - 86 
km SE Mexico City, N slope Mt. Popo¬ 
catepetl, 13500 ft (8 2 6c? TCWC); 36. - 
Paso de Cortes, 3700 m (2 2 lc? CB); 38.- 
Mt. Popocatepetl, 12900 ft * (1 2 USNM). 
Puebla: 37. - 6 km N Paso de Cortes, 3830 
m (9 2 6c? UAMI); 37. - 5 km N Paso de 
Cortes, 3810 m (8 2 7 c? UAMI); 37. - 4 
km N Paso de Cortes, 3710 m (18 2 11c? 
UAMI); 37. - 2.5 km N Paso de Cortes, 
3670 m (22 lc? UAMI); 38. - 15 km NW 
San Martin, Rio Otalti, 8700 ft (3 2 5 c? 
TCWC). 

Thomomys umbrinus martinensis. M = San 
Martin Texmelucan (47). Holotype: Fe¬ 
male, age 2 (55622 USNM).- Puebla: 39. 

— San Luis Coyotcingo (through San Pedro 
Coxtocan Ranch) (1 2 USNM); 40. - Ha¬ 
cienda San Pedro Coxtocan, km 96.5 Pue.- 
Mex. Hgwy. (18 2 26c? USNM); 41. - San 
Martin (Texmelucan], 7400 ft * (2 2 
USNM). 

Thomomys umbrinus albigularis. A = Sierra 
of Pachuca, El Chico (81). Holotype: Fe¬ 
male, age 2 (5188 USNM).- Hidalgo: 42. 

— 1 km E Hacienda Tepozan, Municipio de 
Almoloya, 2700 m (lc? UAMI); 43. - 8 
km E Singuilucan, 2760 m (lc? UAMI); 
44. - Mirasol, 2530 m (3 2 3 c? UAMI); 45. 

- Cuatro Palos, 2500 m (35 2 18c? UAMI); 
46. - Real del Monte, 9500 ft (2 2 lc? 
USNM); 47.- Parque Nacional El Chico, 
2860 m (1 2 UAMI); 47. - El Chico, Sierra 
de Pachuca, 9800 ft * (3 2 USNM); 48. - 
Tulancingo, Sierra de Pachuca, 8500 ft (1 2 
USNM);^49. - 5 mi E Tulancingo, 7400 ft 
(6 2 5 c? KU); 51. - 4 km NW El Carmen, 
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Municipio Almoloya, 2680 m (1 2 UAMI). 
Puebla: 50. — 10 km E, 1 km S Ahuazo- 
tepec, 2300 m (1 2 UAMI). 

Thomomys umbrinus ssp. X = Sierra de 
Tlaxco (369). Tlaxcala: 52. — Acopinalco 
del Penon, 2730 m (3 2 2c? UAMI); 52. - 
Acopinalco del Penon, 2800 m (1 6 
UAMI); 53. - Ejido el Penon, 2960 m 
(17 2 21c? UAMI); 53. - 1 km E Ejido el 
Penon, 2950 m (8 2 1c? UAMI); 54. - El 
Final de la Senda, 2700 m (15 2 3 c? 
UAMI); 54. - El Final de la Senda, 2750 
m (23 2 16c? UAMI); 54. - 1 km S El 
Final de la Senda, 2700 m (12 1c? UAMI); 
54. — Final de la Senda, lfmite Puebla- 
Tlaxcala, 15 km N Tlaxco, 2856m (12 3 c? 
UAMI); 54. - Lfmite Puebla-Tlaxcala, 15 
km N Tlaxco, 2750 m (6 2 5 c? UAMI); 54. 

- Lfmite Puebla-Tlaxcala, 15 km N Tlax¬ 
co, 2750 m (6 2 2c? UAMI); 54. - Lfmite 
Puebla-Tlaxcala, 15 km N Tlaxco, 2856 m 
(7 2 1c? UAMI); 54. - 15 km N, 3 km E 
Tlaxco, 2865 m (12 UAMI); 54. - 2 km 
N, 1 km E Tlaxco, 2880 m (5 2 lc? 
UAMI); 54. - 3 km N, 1 km E Tlaxco, 
2640 m (12 5 c? UAMI); 54. - 3 km N, 1 
km E Tlaxco, 2820 m (87 2 51c? UAMI); 
54—5 km N, 3 km E Tlaxco, 2600 m (1 2 
1 6 UAMI); 54. - 5 km N, 3 km E Tlaxco, 
2920 m (1 2 1<? UAMI); 54. - 5 km N, 3 
km E Tlaxco, 2960 m (2 2 lc? UAMI); 54. 

- 6 km N, 2 km E Tlaxco, 2770 m (2 2 
26 UAMI); 54. - 8 km N Tlaxco, 2820 m 
(1 2 UAMI); 56. - Paso Ancho, 12.5 km 
N, 6 km E Tlaxco, 2620 m (3 2 2 c? 
UAMI); 57. - Villarreal, 3000 m (26 2 
14c? UAMI). Puebla: 55. - Cienega Larga, 
12 km S, 5 km W Chignahupan, 2640 m 
(15 2 12c? UAMI); 55. - Agua Santa, 13.5 
km S, 5 km W Chignahuapan, 2740 m (1 2 
UAMI). 

Thomomys umbrinus ssp. I = Monte La Mal- 
inche (90). Tlaxcala: 58. — La Junta (lc? 
USNM); 59. - Cd. Industrial Xicotencatl, 
2510 m (12 1c? UAMI); 60. - 1 km S Sta. 
Cruz Techachalco, 2080 m (12 UAMI); 
61. — 7 km S, 11 km E Apizaco, 2600 m 
(2 2 UAMI); 61. - 8 km S, 8 km E Api¬ 
zaco, 2600 m (6 2 4c? UAMI); 61. - 8.5 
km S, 9.5 km E Apizaco, 2620 m (7 2 lc? 
UAMI); 61. - 8 km S, 11.5 km E Apizaco, 
2600 m (5 2 8 c? UAMI); 62. - Monte La 
Malinche, 2.5 km S, 11 km W Huamantla, 


3000 m (38 2 13 c? UAMI); 63. - S slope 
Malinche volcano, 3800 m (1 2 USNM); 
64. - San Felipe Tenextepec (12 USNM). 
Thomomys umbrinus orizabae. R = Pico de 
Orizaba Volcanoe (72). Holotype: Female, 
age 3 (53616 USNM).- Puebla: 65. - Mi¬ 
guel Hidalgo y Costilla, 3260 m (24 2 20 c? 
UAMI); 66. - 16 km NE San Andres, W 
slope Orizaba, 11,000 ft (12 TCWC); 66. 
— 12 km NNE San Andres, NW slope Mt. 
Orizaba, 10,000 ft (2 2 TCWC); 66. - 12 
km NNE San Andres, W slope Orizaba, 
10,000 ft (2 2 TCWC); 67. - 6.5 km S, 
10.5 km E Tlalchichuca 3100 m (lc? CB); 
67. — 6.5 km S, 10 km E Tlalchichuca, 
3100 m (12 1c? CB); 68. - 10 km N, 12.4, 
E Cd. Serdan, 3340 m (3 2 lc? CB); 69. - 
Mt. Orizaba, 9500 ft * (102 5<? USNM); 
70. — 2.5 km N, 15.5 km E Ciudad Serdan, 
3280 m (1 2 CB). 

Thomomys umbrinus umbrinus. S = Boca de 
Monte (9). Veracruz: 71. — Boca del Mon¬ 
te, 7800 ft * (4 2 4c? USNM); 72. - 2 km 
N, Los Jacales, 4500 ft (1 2 KU) 
Thomomys umbrinus ssp. T = Totoltepec (5). 
Mexico: 73. - 4.6 km SE Totoltepec, 2540 
m (2 2 3 c? UAMI). 

Thomomys umbrinus ssp. D = La Soledad 
(12). Tlaxcala: 74. — 3 km S La Soledad 
2900 m (102 2 c? UAMI). 

Thomomys umbrinus ssp. Y = Oyameles (5). 
Puebla: 75. - 10 Km E Oyameles, 2500 m 
(12 26 UAMI); 75. - 7 km E Oyameles, 
2550 m (12 3c? UAMI); 76. - 4.2 km W 
Oyameles (12 CB). 

RESULTS 

Nongeographic Variation 

Individual Morphometric Variation 

In general, CVs of the external measure¬ 
ments averaged greater for all ages in both 
sexes, followed by the mandible measure¬ 
ments and then by those of the braincase and 
rostrum (table 2). CVs ranged between 1.9 
for LBC in old females (24) and 14.5 for 
LVT in young males (c?2). Individual varia¬ 
tion differed according to sex with females 
showing less variation than males, except for 
the oldest females (24) whose variation was 
slightly greater than in males of the same age 
(c?4). Indeed, old females (24) exhibited a 
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TABLE 2 

Average Coefficients of Variation (CVs) in a 
Sample of 356 Specimens of Thomomys umbrinus 
from Sierra de Tlaxco in the State of Tlaxcala 

CVs were calculated in 4 external, 22 cranial, and 11 
jaw measurements by sex and age category (juvenile = 
1; subadult = 2; adult = 3; old = 4). Sample sizes (n) are 
shown for each age group; mean and total are given for 
each sex. 


Age 

External 

Cranial 

Jaw 

n 



Females 



1 

7.4 

4.4 

5.0 

22 

2 

8.1 

5.0 

5.7 

104 

3 

8.2 

4.2 

4.9 

80 

4 

7.7 

5.3 

6.7 

8 

Mean 

(Total) 

7.8 

4.7 

5.6 

(214) 



Males 



1 

10.1 

5.4 

8.7 

19 

2 

8.4 

5.0 

7.1 

42 

3 

7.6 

4.1 

4.0 

57 

4 

8.3 

4.1 

6.5 

24 

Mean 

(Total) 

8.6 

4.6 

6.6 

(142) 


greater individual variation in all of their 
measurements than females of the other age 
groups (22—3). Among the latter, subadult 
females (22), revealed a greater overall var¬ 
iation than the young (21) and the adult 
(23) females. Males differed in that the 
young ( 6 \ ) were more variable in all mea¬ 
surements, subadult (c?2) and old (64) males 
showed an overall comparable variation, and 
adult (63) males showed a smaller variation. 

In external measurements (table 2), young 
females (21) showed less variation, followed 
by old females (24), subadults (22), and 
adults (23). In skull characters, adult females 
(23) showed less variation followed by 
young (2 1), subadult (22), and old (24) fe¬ 
males. Adult males (63) showed less varia¬ 
tion in external and mandibular measure¬ 
ments than old (64), subadult (62), and 
young (c?l) males, respectively. However, 
they showed the same lower variation in 
measurements of the braincase as the old 
males (64) as compared to subadult (62) 
and young ( 6 1) males, respectively. 

In both sexes and in all ages, the most var¬ 
iable external characters were LVT (CV = 


6.1-4.6) and LOE (CV = 8.0-12.2), while 
HFT (CV = 5.7-8.0) and TOL (CV = 3.7- 

8.8) were the least variable. For the skull, 
BOP (CV = 6.8-9.1), LMX (CV = 5.0-6.3), 
and LNA (CV = 5.3-6.9) reach a greater 
variation between the females of all ages. 
LNA (CV = 7.4,7.5) together with WRO 
(CV = 7.7, 5.9) and HMC (CV = 7.4, 5.8) 
were also the most variable characters in 
both young (<J1) and subadult males (62). 
Whereas in females, the skull measurements 
with smaller variation were LBC (CV = 1.9— 

2.9) , TLS (CV = 2.0—3.2), and SBP (CV = 

2.5— 4.1), in males it was BAB (CV = 2.5— 
3.7). For the jaw, measurements with greater 
variation in females were AWD (CV = 7.6— 
9.2), DHI (CV = 5.2-8.6), and BJT (CV = 

5.5— 7.4), while in males it was also AWD 
(CV = 7.0-12.1), besides AHM (CV = 9.4- 
10.0), and BML (CV = 7.7-11.9). Mensu¬ 
rable characters with smaller CVs were BCP 
(CV = 2.8—7.7) and GLD (CV = 2.9-6.2) 
in both sexes; although, young males (<?1) 
showed a greater variation. 

Ontogenetic Variation 

Table 3 shows the ontogenetic variation 
among subadult (age 2), adult (age 3), and 
old specimens (age 4) according to a single 
ANOVA. Results show that the age catego¬ 
ries separated in a high proportion of vari¬ 
ables. Differences in patterns are also shown 
by the sexes: in males, fewer variables 
showed no differences between ages (LVT, 
LOE, and BAM) than in females (HFT, LVT, 
ICO, LPF, BAM, BOP, HBP, and BML). Es¬ 
pecially, old females (age 4) segregated more 
from other females than males of the same 
age with respect to other males (table 3). 

In the analyses in which ages 1—4 were 
considered collectively, variables with no 
differences among females were LVT, HFT, 
and ICO; only ICO and BOP showed no dif¬ 
ferences among males of different ages. 
Among younger specimens (age 1), ICO dif¬ 
fered more than that of the other ages. Young 
specimens segregated in almost every vari¬ 
able, regardless of their sex, and, therefore, 
they were excluded from all of the later anal¬ 
yses. In general, in both analyses, variables 
showing discontinuity between age groups, 
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TABLE 3 

Effect of Age on 37 Variables (Var) in a Sample of 306 Specimens of Thomomys umbrinus 
from Sierra deTlaxco in the State of Tlaxcala, According to an Analysis of Variance (F) 

Within Duncan's test, a hyphen indicates no differences between the means of the involved age groups (subadults 
= 2; adults = 3; old = 4), and a comma indicates differences. All tests were significant (0.0001 =s p =s 0.05) unless 
marked by with asterisk. For definitions of variables, see the Measurements section in the Introduction. 


Var 

Males 

Females 

Var 

Males 

Females 

F 

Duncan 

F 

Duncan 

F 

Duncan 

F 

Duncan 

TOL 

33.0 

4, 3, 2 

19.3 

4, 3,2 

LVT 

3.6 

4-3, 3-2 

*1.6 

4-3-2 

HFT 

8.0 

4-3, 2 

4.1 

4-3, 3-2 

LOE 

3.5 

4-3, 3-2 

4.3 

4-3, 2 

GLS 

82.1 

4, 3, 2 

76.9 

4, 3, 2 

LNA 

46.0 

4, 3, 2 

55.4 

4, 3, 2 

LMX 

55.6 

4, 3, 2 

40.1 

4, 3, 2 

LFO 

12.7 

4-3, 2 

9.6 

4, 3-2 

LZA 

67.5 

4, 3, 2 

54.0 

4, 3, 2 

LBC 

70.5 

4, 3, 2 

53.9 

4, 3, 2 

WRO 

48.0 

4, 3, 2 

41.7 

4, 3, 2 

ICO 

5.9 

4, 3-2 

*1.6 

4-3-2 

AZB 

100.8 

4, 3, 2 

56.8 

4, 3, 2 

MZB 

97.7 

4, 3, 2 

49.5 

4, 3, 2 

PZB 

114.1 

4, 3, 2 

49.7 

4, 3, 2 

SBP 

8.5 

4-3, 2 

14.7 

4, 3, 2 

LDI 

105.6 

4, 3, 2 

66.1 

4, 3, 2 

LMT 

7.4 

4, 3-2 

13.3 

4, 3, 2 

LPF 

31.6 

4-3. 2 

6.7 

4-3, 3-2 

BAM 

8.3 

4-3, 2 

4.6 

4-3, 3-2 

BOP 

5.2 

4-3, 2 

3.1 

4-3, 3-2 

BAB 

67.6 

4, 3,2 

39.4 

4, 3, 2 

HRO 

38.5 

4, 3, 2 

26.2 

4, 3, 2 

HMC 

101.5 

4, 3, 2 

79.1 

4, 3, 2 

MHS 

57.0 

4, 3, 2 

62.7 

4, 3, 2 

HBP 

7.6 

4-3, 2 

5.1 

4-3, 3-2 

BCP 

7.1 

4-3, 2 

18.3 

4-3, 2 

BJT 

*2.7 

4-3, 3-2 

14.0 

4, 3-2 

AWD 

23.3 

4, 3, 2 

19.2 

4, 3, 2 

GWM 

18.9 

4, 3, 2 

68.6 

4, 3, 2 

MWD 

6.9 

4-3, 2 

28.2 

4, 3, 2 

GLD 

69.5 

4, 3, 2 

68.0 

4, 3, 2 

LCP 

44.7 

4, 3, 2 

49.0 

4, 3, 2 

BML 

10.0 

4-3, 2 

*1.5 

4-3-2 

DHI 

31.8 

4, 3, 2 

41.0 

4, 3, 2 

AHM 

8.2 

4, 3, 2 

52.3 

4, 3, 2 

PHM 

69.9 

4, 3, 2 

53.1 

4, 3, 2 







are related to body length and to the length 
and width of the rostrum and mandible. 

Secondary Sexual Dimorphism 

In Thomomys umbrinus (table 4), signifi¬ 
cant sexual dimorphism (p < 0.05) was more 
evident among adult (83.8% of the variables) 
and old specimens (73% of the variables) 
than among the subadults (27% of the vari¬ 
ables). Except for ICO, where subadult fe¬ 
males surpassed males of the same group, 
males exceeded females in all other charac¬ 
ters with significant sexual dimorphism. 

Males were also larger than females in 
several variables with no significant sexual 
dimorphism (old specimens = 35.1% of var¬ 
iables; adult specimens = 13.5%, and sub¬ 
adult specimens = 27%). There were only 
six measurements where females were sig¬ 
nificantly larger than males of the corre¬ 
sponding age group: subadult females were 
larger in BOP, MWD, BAM, BJT, and AHM; 
adult females in BOP; and old females in 
BML. 


Combined Effect of Age and Sexual 
Dimorphism 

Analyses of the combined effect of age 
and sexual variation supported previous re¬ 
sults (table 5). Age groups were separated in 
every jaw (11 variables), 20 cranial, and two 
external measurements, allowing the differ¬ 
entiation among the three considered ages. 
Exceptions were LVT, with no differences 
between adult and subadult specimens; LOE 
with no differences between old and adult 
specimens; ICO with nonsignificant differ¬ 
ences among the age groups, although sub¬ 
adults were broader than the other groups; 
and BOP which showed no differences be¬ 
tween old and adult individuals and between 
adult and subadults. 

Sexual dimorphism was also present in al¬ 
most every variable (35 out of 37), except 
for ICO and BOP where females were non- 
significantly larger than males (table 5). In¬ 
teraction between age and sexual dimor¬ 
phism was positive in one-third of all char¬ 
acters (32%), including one body measure- 
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TABLE 4 

Sexual Dimorphism in 306 Specimens of Thomomys umbrinus 
from Sierra de Tlaxco in the State of Tlaxcala 

The ANOVA was performed on 37 variables (Var) using sex and age groups. F values (F) and Tukey’s tests (Tu) 
are shown by age group (subadult = 2, n = 141; adult = 3, n = 133; old = 4, n = 32). Within the latter, > indicates 
that males were significantly (0.0001 =s p =s 0.05) larger than females; a indicates that males were nonsignificant- 
ly larger than females; < indicates that females were significantly larger than males; and v indicates that females 
were nonsignificantly larger than males. 


Var 



Age group 



Var 



Age group 



2 


3 


4 


2 


3 


4 


F 

Tu 

F 

Tu 

F 

Tu 

F 

Tu 

F 

Tu 

F 

Tu 

TOL 

2.7 

A 

33.3 

> 

14.0 

> 

LVT 

0.7 

A 

6.9 

> 

1.7 

A 

HFT 

3.0 

A 

36.1 

> 

12.0 

> 

LOE 

11.9 

> 

10.7 

> 

3.3 

A 

GLS 

7.3 

> 

120.6 

> 

17.0 

> 

LNA 

3.4 

A 

44.2 

> 

16.4 

> 

LMX 

4.6 

> 

52.2 

> 

16.2 

> 

LFO 

0.9 

A 

14.7 

> 

3.0 

A 

LZA 

1.0 

A 

57.0 

> 

2.9 

A 

LBC 

6.5 

> 

96.2 

> 

7.9 

> 

WRO 

3.6 

A 

51.1 

> 

15.1 

> 

ICO 

4.3 

< 

6.2 

> 

3.0 

A 

AZB 

17.7 

> 

167.6 

> 

22.1 

> 

MZB 

21.5 

> 

147.5 

> 

17.8 

> 

PZB 

21.9 

> 

141.8 

> 

12.8 

> 

SBP 

0.0 

A 

3.0 

A 

2.0 

A 

LDI 

9.4 

> 

169.2 

> 

9.0 

> 

LMT 

0.1 

A 

1.6 

A 

4.1 

> 

LPF 

3.2 

A 

25.8 

> 

0.5 

A 

BAM 

0.0 

V 

2.5 

A 

1.0 

A 

BOP 

0.8 

V 

0.3 

V 

1.6 

A 

BAB 

6.5 

> 

86.5 

> 

16.0 

> 

HRO 

1.0 

A 

15.4 

> 

3.0 

A 

HMC 

8.3 

> 

89.7 

> 

9.5 

> 

MHS 

1.3 

A 

40.2 

> 

15.3 

> 

HBP 

2.9 

A 

12.5 

> 

4.5 

> 

BCP 

3.8 

A 

10.4 

> 

12.7 

> 

BJT 

1.1 

V 

0.0 

A 

4.9 

> 

AWD 

2.8 

A 

26.4 

> 

10.3 

> 

GWM 

1.9 

A 

18.1 

> 

19.6 

> 

MWD 

0.1 

V 

2.7 

A 

11.0 

> 

GLD 

3.1 

A 

69.4 

> 

5.2 

> 

LCP 

2.4 

A 

37.0 

> 

9.7 

> 

BML 

0.6 

A 

11.9 

> 

0.1 

V 

DHI 

0.7 

A 

49.6 

> 

5.0 

> 

AMH 

0.0 

V 

4.5 

> 

8.8 

> 

PHM 

2.8 

A 

69.0 

> 

8.4 

> 









ment (LOE), nine cranial measurements 
(TLS, LZA, LBC, ICO, PZB, LDI, LPF, 
BAB, and HMC), and three jaw measure¬ 
ments (GLD, DHI, PHM). 

Finally, in the analysis of the variance 
components where only age, sexual dimor¬ 
phism, and their interaction were considered, 
almost half of variation was due to these fac¬ 
tors (46.91%). Therefore, 25.94% of varia¬ 
tion is due to age effect, 19.42% to second¬ 
ary sexual dimorphism, and 1.55% to the in¬ 
teraction between them. The rest of variation, 
not specified as a variable and included with¬ 
in the standard error (53.09%), might be due 
to other factors such as fossorial niche, type 
of soil, food habits, etc. 

Geographic Variation 
Analyses of Variance 

F statistics and Duncan’s Test for Multiple 
Means are available upon request. In general. 


results of these analyses demonstrated the 
absence of an overall pattern of variation re¬ 
lated to the geographic distribution of the 
samples in all age classes of both sexes along 
the MTB (table 6). On the contrary, with the 
exception of Thomomys u. pullus, the general 
pattern seemed to be that of morphometric 
similarity among the studied populations re¬ 
gardless of their origin (table 6). 

Indeed, every analysis was significant (p 
< 0.05) for the majority of the 24 variables 
examined and F ranged between 0.94 (HRO) 
and 14.58 (DHI), both for the male sub¬ 
adults. Duncan’s tests showed a highly het¬ 
erogeneous series among the group locality 
(GF) subgroups that are formed and varied 
from one sex to the other, from one age to 
another and within a single age and sex. 

Although the results pointed out differenc¬ 
es between the extreme GFs (with greater 
and lesser magnitude in one variable), in 
general, the intermediate GFs formed a con- 
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TABLE 5 

Sexual Dimorphism (S), Effect of Age (A), and Interaction Between These Two Factors (SA) 
in 306 Specimens of Thomomys umbrinus from Sierra de Tlaxco in the State of Tlaxcaia 

Analyses were performed on 37 variables (Var) obtained from 141 subadult (age group 2), 133 adult (age group 3), 
and 32 old (age group 4) specimens. F values (F) and Duncan’s tests for Multiple Means (Duncan) are shown. 
Unless marked by an asterisk, analyses were significant (0.0001 =s p =s 0.05). To = F value for the whole analysis. 
See tables 3 and 4 for notation. 


Var 


F 



Duncan 

Var 


F 



Duncan 

To 

S 

A 

SA 

S 

A 

To 

S 

A 

SA 

S 

A 

TOL 

37.7 

94.8 

46.3 

*0.5 

> 

4, 3,2 

LVT 

5.4 

16.2 

4.9 

*0.5 

> 

4, 3-2 

HFT 

19.1 

70.5 

11.1 

*1.4 

> 

4, 3, 2 

LOE 

8.6 

27.2 

3.7 

4.2 

> 

4-3, 2 

GLS 

11.3 

244.5 

161.3 

4.7 

> 

4, 3, 2 

LNA 

68.9 

138.3 

102.0 

*1.1 

> 

4, 3, 2 

LMX 

66.9 

145.2 

93.7 

*0.9 

> 

4, 3, 2 

LFO 

16.5 

36.1 

22.4 

*0.7 

> 

4, 3, 2 

LZA 

73.5 

114.3 

119.6 

7.0 

> 

4, 3, 2 

LBC 

87.2 

176.6 

123.1 

6.5 

> 

4, 3, 2 

WRO 

65.8 

144.0 

89.6 

*2.9 

> 

4, 3, 2 

ICO 

3.4 

*3.2 

*0.8 

6.1 

V 

2-3-4 

AZB 

110.0 

253.1 

145.8 

*2.6 

> 

4, 3, 2 

MZB 

101.5 

230.2 

136.0 

*2.7 

> 

4, 3, 2 

PZB 

105.3 

228.0 

144.9 

4.5 

> 

4, 3, 2 

SBP 

12.5 

16.1 

22.9 

*0.2 

> 

4, 3, 2 

LDI 

118.0 

236.2 

168.6 

8.3 

> 

4, 3,2 

LMT 

11.6 

17.2 

20.3 

*0.3 

> 

4, 3, 2 

LPF 

22.5 

46.2 

28.6 

4.5 

> 

4, 3, 2 

BAM 

6.2 

9.5 

10.4 

*0.2 

> 

4, 3, 2 

BOP 

3.0 

*0.2 

7.1 

*0.3 

V 

4-3, 3-2 

BAB 

80.2 

191.4 

101.1 

3.7 

> 

4, 3, 2 

HRO 

35.9 

52.2 

62.5 

*1.1 

> 

4, 3, 2 

HMC 

114.2 

198.8 

180.9 

5.1 

> 

4, 3, 2 

MHS 

75.0 

134.5 

119.7 

*0.6 

> 

4, 3, 2 

HBP 

12.3 

34.9 

12.9 

*0.4 

> 

4, 3, 2 

BCP 

18.3 

46.5 

23.5 

*0.6 

> 

4, 3, 2 

BIT 

8.1 

6.8 

14.3 

*2.5 

> 

4, 3,2 

AWD 

32.7 

77.1 

42.4 

*0.6 

> 

4, 3,2 

GWM 

45.6 

84.7 

71.7 

*0.0 

> 

4, 3, 2 

MWD 

18.3 

26.4 

31.1 

*1.5 

> 

4, 3,2 

GLD 

85.8 

140.3 

140.1 

4.4 

> 

4, 3, 2 

LCP 

60.2 

108.4 

95.2 

*1.1 

> 

4, 3, 2 

BML 

5.9 

9.2 

7.6 

*2.5 

> 

4, 3, 2 

DHI 

49.2 

95.1 

70.9 

4.5 

> 

4, 3, 2 

AMH 

22.0 

30.9 

38.8 

*0.9 

> 

4, 3, 2 

PHM 

78.6 

143.5 

121.1 

3.6 

> 

4, 3, 2 









tinuum between the extremes. In fact, there 
were several measurements in which there 
were no differences among GLs: LPF and 
BAM in adult females; AZB and MHS in old 
females; LNA, LZA, and LBC in subadult 
males; TLS, LNA, LZA, ICO, LPF, MHS, 
and AHM in adult males. 

In contrast to the unclear general pattern, 
some individual groups from some sampled 
localities consistently segregated from those 
in other GLs for several variables. In GL P 
(T. u. pullus ), subadult females segregated 
from other GLs for its small size in TLS, 
AZB, BAM, BAB, MHS, HMC and subadult 
males separated for their smaller size in 
BAM and BAB. In GL I (T. umbrinus ssp. 
from La Malinche Volcano), subadult fe¬ 
males segregated according to larger size in 
AHM and old females for larger size in LBC, 
HRO, GLD, and DHI. In GL V (T. u. vul- 
canius ), adult females segregated according 
to larger size in SBP, old females in DHI 
(following GL MM), and old males in BAB. 
In GL R (T. u. orizabae) adult males were 


smaller-sized in LD1. In GL M ( T. u. marti- 
nensis) old males were larger-sized in LCP. 
Finally, in GL N (71 umbrinus ssp. of the 
Sierra de San Andres) old males were larger- 
sized in LZA, HMC, and GLD. 

Principal Components Analysis and 
Minimum Spanning Tree 

Except for the GL P representing 7. u. pul¬ 
lus (fig. 3), all GLs overlapped broadly in 
multispace (figs. 3, 4). Some of the geo¬ 
graphically nearest GLs tended to group to¬ 
gether (i.e., GLs A and X); however, this is 
not a rule in every analysis (i.e., minimum 
spanning trees, fig. 4). Moreover, no consis¬ 
tent geographic pattern related to origin of 
samples was observed among the GLs 
groupings according to age or sex (figs. 3, 
4). 

Distribution patterns of the GLs in multis¬ 
pace were similar between the PCA carried 
out on the individuals (fig. 3) and the PCA 
carried out on the sample means of the GLs 
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TABLE 6 

Studied Populations of Thomomys umbrinus Along the Mexican Transvolcanic Belt Grouped 
by Age Category According to Selected Variables in a Multiple Analysis of Variance 

Mean size of variables decreases from the top to the bottom. In the Duncan analysis, vertical lines indicate absence 
of significant differences (p «= 0.05) among populations and horizontal lines indicate significant segregation from 
the other populations. See table 1 for locality names. 

Age group 2 _ Age group 3 _ _ Age group 4 _ 

Females Males Females Males Females Males 


O 

V 

E 

I 

B 

R 

A 

L 

N 

C 

D 

X 

P 


V 

I 

E 

C 

B 

M 

X 

R 

A 

P 

F= 1.42 


Greatest Length of Skull (GLS) 


O 

R 

V 

L 

M 

A 

C 

X 


o 

R 

A 

B 

X 

C 

F = 1.48 


V I 

I 

o 

R 

X 

A l 

F = 8.59 


V I 
R I 
O 
I 

A 

M I 

X ! 

N 

F = 12.15 


I 

B 

E 

A 

O 

R 

D 

X 

N 

L 

C 

V 

P 

F = 5.63 


Length to Coronoid Process (LCP) 


I 


B 

R 




C 


I 




X 


M 




A 


A 



A 

R 


E 


O 



B 

I 


M 


L 


R 

X 


V 


X 


X 

o 


P 


c 


O 

V 

R 

V 

C 

A 


F = 

= 1.52 

F = 

= 3.88 

F = 2.39 

F = 


M- 

R 

I 

A 

V 

O 

X 

N 

F = 6.94 


I 

O 

D 

A 

B 

X 

E 

V 

N 

C 

L 

R 

P- 

F = 9.86 


Greatest Length of Dentary (GLD) 


I 

A 

X 

M 

C 

E 

B 

R 

P 

V 

F = 2.23 


T 

O 

M 

A 

R 

X 

C 

L 

V 

F = 4.73 


A 

X 

B 

O 

R 

C 

F = 2.39 


I - 

V 

X 

o 

R 

A 

F = 7.64 


I 

V 

R 

A 

M 

X 

O 

N- 

F = 5.94 
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(Continued) 


Age group 2 


Females 


Males 


Age group 3 


Age group 4 


Females 


Males 


Females 


Males 


D 

1 

O 

A 

B 

V 

E 

X 

R 

C 

N 

L 

P 


F = 12.81 


F = 2. 30 


Anterior Zygomatic Breadth (AZB) 


X 


O 



E 


I 



B 


R 




M 


A 


A 

I I 

V 


V 


X 


V 

A 


C 


B 


o 

C 


M 


O 


X 

R 


X 

C 

A 

P 

L 

R 

1 R 


F = 3.56 


F = 3.28 


F= 1.91 


V 

M 

A 

I 

O 

R 

X 

N 


F = 2.16 


Length of Braincase (LBC) 


I 

B 

O 

E 

D 

L 

V 

X 

A 

R 

C 

N 

P 


F = 5.13 


E 

B 

I 


I 

R 


V 

O 




R 


M 

L 


A 

I - 

V 


C 

V 


B 

V 

0 


X 

M 


O 

o 

A | 

R 

C 


X 

R 


X 

A 

A 


R 

X 


M 

P 

X 

C 

A 

N 


F = 1.36 

F = 

= 8.15 

F = 1.07 

F = 

= 7.0 

F = 


(fig. 4). However, both analyses differed (ta¬ 
ble 7) in the most important variables for PC 
1—111, the percentage of variation explained 
by PC 1-111, and the relation of PC-1 to either 
size or shape. 

In the PCA on individuals (table 7), 65.7% 
of variation was explained by the first three 
PCs. The highest amount of variation ex¬ 
plained was reached in the analysis for sub¬ 
adult males (72.3%) and the lower in both 
adult males and females, 58% and 58.8%, 
respectively (fig. 3). PC 1 explained more 
than 50% of the total variation, except in 
both adult specimens (males, 40%; females, 
42.7%). P 11 explained from 7.2% in subadult 
males to 11.1% in old females, and PC 111 


explained from 4.9% in old males to 8.5% in 
adult males. 

All coefficients had positive sign in PC 1, 
regardless of sex or age (fig. 3, table 7); 
therefore, in such analyses PC 1 constituted 
a component related to the size of the indi¬ 
viduals (Pimentel, 1979). Coefficients for PC 
11 and PC 111 had different sign in all PC As; 
therefore these components are related to 
shape (Pimentel, 1979). According to their 
coefficients, the first five most important cra¬ 
nial measurements for PC 1 were TSL, LCP, 
GDL, LBC, and AZB; for PC 11 were LCP, 
LFO, BAB, TSL, and GDL; and for PC 111 
were LCP, LFO, BAB, LNA, and HBP. 

In the PCA on the GL means (fig. 4, table 
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PCI 


Fig. 3. Bidimensional distribution of subadult 
females in the principal components (PC) III and 
I. Thomomys umbrinus pullus (P) is separated 
from other examined populations along the Mex¬ 
ican Transvolcanic Belt by both its small size and 
cranial morphology. Letters indicate the centroids 
of each examined sample as mentioned on table 
1. The analysis was made on the individuals. 


7), 79% of variation was explained by the 
first three PCs. In this analysis (fig. 4), the 
highest amount of variation was explained in 
the analysis for adult males (89.4%) and the 
lowest in that for females of the same age 
(69.5%). PC I explained from 36.9% (adult 
females) to 74.1% (adult males) of the total 
variation among samples; PC II from 9.4% 
(adult males) to 33.9% (subadult females); 
and PC III from 4.9% (old males) to 15.4% 
(old females). 

Since all coefficients had positive signs in 
PC I for old males (fig. 4, table 7), this func¬ 
tion represents size tendencies among GLs 
(Pimentel, 1979). Coefficients in the other 
analyses had different signs for all PCs (fig. 
4, table 7); therefore, these components are 


related to shape (Pimentel, 1979). In PC I, 
the most important measurements for old 
males were TLS, GDL, LCP, LBC, and AZB, 
while for the other specimens they were 
TLS, LCP, PHM, GDL, AZB. In PC II the 
most important measurements were TLS, 
LCP, LNA, GDL; LFO, and BAB, and in PC 
Ill LFO, LMX, PHM, GDL, and LCP. 

Comparison of patterns in plots of figure 
4 revealed no consistent pattern among the 
outstanding variables in relation to sex, age 
group, or GL, because they were not always 
the same or in the same sequence (table 7). 
However, it is evident that the coefficients of 
the five most important variables surpassed 
the others considerably in every case and that 
there were variables with almost no impor¬ 
tance at all (table 7). 

The majority of the coefficients in the PCA 
on individuals, as in the PCA on the means, 
showed a positive sign that indicates an in¬ 
crease in the dimensions of these variables 
among the studied populations (table 7). On 
the other hand, the variables that showed a 
negative coefficient can be interpreted as hav¬ 
ing a decrease in magnitude among the pop¬ 
ulations (Pimentel, 1979). Such variables 
were LFO in subadult females; BAB and 
HBP in adult females; AZB and SBP in old 
females; LNA, LMX, LFO, WRO, and LPF 
in subadult males; and TLS, LBC, WRO, 
AZB, LPF, BAM, BAB, HRO, HMC, GLD, 
and LCP in adult males (table 7). 

It can be pointed out that all of the vari¬ 
ables of the first analysis were included in 
those outstanding in the second analysis, ex¬ 
cept for SBP, LZA, and HBP (table 7). Re¬ 
garding size, these variables (especially TLS, 
LFO, LCP, and GLD) allowed the distinction 
of the T. u. pullus individuals from the rest 
due to their shorter dimensions. In shape, this 
subspecies was also distinguished from other 
examined T. umbrinus by a lesser expansion 
of the zygomatic arches (AZB, MZB, PZB) 
and by a general smoother appearance of the 
skull. 

According to the minimum spanning tree 
(MST) projected onto the PCA (fig. 4), some 
populations showed morphophenetic rela¬ 
tions with more than one population (GLs X, 
O, B, A, R, N, M, L, C, I), while others 
linked with only one (P, T, D, E, S, Y). Phe- 
netic distances were shorter between some 
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TABLE 7 

Percentage of Variation Explained by the First Three Principal Components (% PC I—HI) and 
Their Coefficients for 24 Cranial Variables (Var) Measured in Thomomys umbrinus 
from 17 Populations Along the Mexican Transvolcanic Belt 

A first principal components analysis (PCA) was computed using all specimens (Individuals) and a second PCA 
with the population averages (Means). Analyses were developed by separated sexes and age groups. The total per¬ 
centage of variation explained by PC I—III (% PC) is given for each analysis. 


Individuals, Females 

% Total: 

Age group 2 = 

65.7% 

Age group 3 = 

58.8% 

Age group 4 = 

68.8% 

Var 

PC I 

PC II 

PC III 

PCI 

PC II 

PC III 

PCI 

PC II 

PC III 

GLS 

0.47 

-0.45 

0.15 

0.55 

-0.44 

0.17 

0.45 

-0.38 

-0.06 

LNA 

0.18 

-0.29 

-0.27 

0.23 

-0.35 

-0.28 

0.21 

-0.09 

0.02 

LMX 

0.16 

-0.22 

-0.24 

0.20 

-0.22 

-0.15 

0.20 

-0.24 

0.10 

LFO 

0.12 

-0.48 

0.44 

0.16 

-0.17 

0.43 

0.16 

-0.48 

0.54 

LZA 

0.21 

-0.02 

-0.11 

0.19 

0.15 

0.15 

0.25 

0.17 

0.16 

LBC 

0.25 

-0.12 

0.02 

0.28 

0.11 

0.20 

0.26 

0.14 

-0.15 

WRO 

0.09 

-0.06 

0.05 

0.13 

0.02 

-0.01 

0.11 

-0.07 

-0.02 

ICO 

0.04 

0.00 

0.08 

0.03 

0.10 

0.10 

0.05 

-0.13 

0.18 

AZB 

0.33 

0.33 

0.21 

0.25 

0.10 

-0.14 

0.22 

0.10 

0.14 

SBP 

0.06 

0.06 

0.14 

0.08 

0.11 

0.22 

0.15 

-0.11 

0.29 

LDI 

0.23 

-0.03 

-0.15 

0.21 

-0.20 

-0.11 

0.24 

-0.09 

0.05 

LMT 

0.08 

-0.00 

-0.02 

0.09 

0.10 

0.01 

0.10 

-0.03 

0.03 

LPF 

0.07 

0.04 

-0.00 

0.10 

0.04 

0.03 

0.07 

0.03 

0.04 

BAM 

0.09 

0.13 

0.16 

0.03 

0.06 

-0.07 

0.06 

-0.10 

0.01 

BAB 

0.22 

0.38 

0.41 

0.14 

0.31 

0.55 

0.23 

-0.13 

0.33 

HRO 

0.08 

-0.02 

0.03 

0.11 

0.07 

0.00 

0.12 

0.07 

0.07 

HMC 

0.19 

-0.01 

-0.08 

0.21 

0.03 

-0.06 

0.21 

-0.02 

-0.09 

MHS 

0.14 

0.09 

0.12 

0.12 

0.03 

-0.09 

0.14 

0.10 

0.05 

HBP 

0.09 

0.13 

0.15 

0.07 

0.08 

0.11 

0.08 

0.19 

0.07 

GLD 

0.32 

0.13 

-0.00 

0.27 

0.21 

-0.09 

0.33 

0.17 

0.22 

LCP 

0.35 

0.21 

-0.55 

0.31 

0.35 

-0.36 

0.24 

0.42 

-0.54 

DHI 

0.05 

-0.00 

-0.04 

0.08 

0.08 

-0.03 

0.10 

0.01 

0.10 

AHM 

0.12 

0.08 

-0.07 

0.11 

0.16 

-0.09 

0.18 

0.12 

0.12 

PHM 

0.20 

0.20 

-0.10 

0.16 

0.41 

-0.21 

0.21 

0.39 

0.07 

% PC 

51.6 

8.4 

5.7 

42.7 

9.1 

7.0 

50.0 

11.0 

7.8 






Individuals, Males 




% Total: 

Age group 2 = 

72.3% 

Age group 3 = 

58.0% 

Age group 4 = 

70.3% 

Var 

PCI 

PC II 

PC III 

PC I 

PC II 

PC III 

PCI 

PC II 

PC III 

GLS 

0.41 

-0.28 

0.36 

0.43 

0.29 

-0.28 

0.51 

-0.41 

-0.12 

LNA 

0.20 

-0.47 

-0.13 

0.24 

0.06 

0.14 

0.20 

-0.01 

-0.38 

LMX 

0.17 

-0.33 

-0.22 

0.23 

0.06 

0.14 

0.24 

-0.04 

-0.10 

LFO 

0.07 

-0.03 

0.44 

0.10 

0.42 

-0.60 

0.21 

-0.46 

0.34 

LZA 

0.20 

-0.24 

0.04 

0.24 

0.10 

-0.02 

0.20 

0.10 

0.36 

LBC 

0.22 

-0.12 

0.27 

0.23 

0.13 

-0.10 

0.27 

-0.09 

0.24 

WRO 

0.07 

-0.05 

-0.01 

0.13 

-0.09 

-0.07 

0.14 

-0.43 

-0.12 

ICO 

0.02 

0.09 

0.14 

0.02 

0.02 

-0.11 

0.03 

-0.06 

0.05 

AZB 

0.33 

0.20 

-0.02 

0.25 

0.01 

0.31 

0.27 

0.18 

-0.45 

SBP 

0.05 

0.10 

0.20 

0.07 

0.25 

0.11 

0.09 

-0.11 

0.03 

LDI 

0.26 

-0.15 

-0.08 

0.21 

0.21 

0.22 

0.19 

-0.09 

-0.06 

LMT 

0.06 

-0.06 

0.09 

0.09 

-0.12 

-0.14 

0.11 

0.01 

0.04 

LPF 

0.07 

0.11 

0.07 

0.06 

0.06 

-0.00 

0.10 

0.03 

0.14 

BAM 

0.05 

0.24 

0.19 

0.04 

-0.02 

-0.03 

0.04 

-0.04 

-0.09 

BAB 

0.19 

0.39 

0.23 

0.17 

0.08 

0.36 

0.17 

0.02 

-0.26 
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TABLE 7 
(Continued) 







Individuals, Males 




% Total: 

Age group 2 = 

: 72.3% 

Age group 3 = 

: 58.0% 

Age group 4 = 

: 70.3% 

Var 

PCI 

PC II 

PC III 

PCI 

PC II 

PC III 

PC 1 

PC II 

PC III 

HRO 

0.08 

-0.01 

0.10 

0.10 

0.01 

-0.08 

0.11 

-0.02 

0.15 

HMC 

0.19 

-0.11 

-0.01 

0.19 

0.05 

0.11 

0.19 

0.02 

0.04 

MHS 

0.16 

0.10 

0.09 

0.16 

0.00 

0.07 

0.12 

0.07 

0.07 

HBP 

0.11 

0.16 

0.24 

0.11 

0.17 

0.30 

0.04 

-0.03 

0.28 

GLD 

0.36 

0.13 

-0.13 

0.35 

-0.14 

0.03 

0.30 

0.19 

0.11 

LCP 

0.39 

0.35 

-0.49 

0.36 

-0.50 

-0.23 

0.28 

0.50 

0.09 

DHI 

0.08 

-0.02 

-0.02 

0.06 

0.07 

0.07 

0.08 

0.02 

0.04 

AHM 

0.11 

0.08 

-0.05 

0.16 

-0.45 

-0.07 

0.13 

0.13 

0.10 

PHM 

0.21 

0.11 

-0.13 

0.19 

-0.16 

0.17 

0.16 

0.46 

0.18 

% PC 

59.1 

7.2 

6.1 

40.0 

9.5 

8.5 

56.7 

8.7 

4.9 






Means, Females 




% Total: 

Age group 2 = 

82.3% 

Age group 3 = 

69.5% 

Age group 4 = 

75.0% 

Var 

PCI 

PC 11 

PC III 

PCI 

PC II 

PC III 

PCI 

PC II 

PC III 

GLS 

0.37 

0.16 

-0.53 

0.56 

0.02 

-0.13 

0.43 

-0.25 

0.40 

LNA 

0.05 

0.10 

-0.02 

0.31 

-0.40 

-0.12 

0.22 

-0.48 

-0.03 

LMX 

0.03 

-0.00 

0.13 

0.28 

-0.30 

-0.41 

0.32 

-0.45 

-0.05 

LFO 

-0.00 

0.23 

-0.43 

0.24 

0.15 

-0.07 

0.09 

0.07 

0.66 

LZA 

0.11 

0.16 

0.26 

0.19 

0.14 

0.14 

0.21 

0.37 

-0.03 

LBC 

0.18 

0.18 

-0.13 

0.42 

0.11 

0.12 

0.24 

0.30 

-0.03 

WRO 

0.08 

0.06 

-0.06 

0.16 

-0.05 

-0.05 

0.14 

0.22 

0.12 

ICO 

0.06 

0.05 

-0.01 

0.03 

0.11 

-0.07 

0.08 

0.08 

-0.02 

AZB 

0.33 

0.34 

0.24 

0.12 

-0.10 

0.15 

-0.03 

-0.01 

-0.12 

SBP 

0.03 

0.12 

0.19 

0.08 

0.12 

-0.09 

-0.07 

-0.03 

0.13 

LDI 

0.19 

0.07 

-0.22 

0.13 

-0.03 

0.09 

0.29 

-0.04 

-0.10 

LMT 

0.13 

0.06 

-0.16 

0.11 

0.12 

-0.03 

0.04 

-0.07 

0.07 

LPF 

0.09 

-0.01 

0.10 

0.12 

0.16 

0.04 

0.11 

0.08 

0.01 

BAM 

0.07 

0.20 

0.06 

0.03 

0.05 

-0.01 

0.03 

0.00 

0.07 

BAB 

0.24 

0.22 

0.22 

-0.04 

0.62 

-0.62 

0.19 

0.03 

0.08 

HRO 

0.05 

0.08 

0.11 

0.07 

0.01 

-0.07 

0.09 

0.15 

0.09 

HMC 

0.17 

0.07 

-0.20 

0.18 

0.18 

0.19 

0.18 

-0.03 

0.12 

MHS 

0.14 

0.14 

-0.03 

0.08 

0.16 

0.22 

0.03 

-0.03 

-0.24 

HBP 

0.09 

-0.02 

0.06 

-0.02 

0.17 

0.16 

0.17 

0.23 

0.04 

GLD 

0.27 

0.23 

0.18 

0.09 

0.34 

0.25 

0.27 

-0.04 

-0.28 

LCP 

0.58 

-0.72 

-0.01 

0.23 

0.09 

0.24 

0.31 

0.22 

-0.16 

DHI 

0.04 

0.02 

0.07 

0.14 

-0.06 

-0.02 

0.12 

0.16 

0.07 

AHM 

0.13 

0.08 

0.13 

0.07 

0.07 

0.01 

0.20 

0.17 

-0.03 

PHM 

0.27 

-0.07 

0.27 

0.03 

0.07 

0.28 

0.28 

-0.01 

-0.37 

% PC 

40.1 

33.9 

8.3 

36.9 

18.0 

14.6 

40.3 

19.3 

15.4 






Means, Males 




% Total: 

Age group 2 = 

74.2% 

Age group 3 = ; 

89.4% 

Age group 4 = ! 

85.7% 

Var 

PC I 

PC II 

PC III 

PC I 

PC II 

PC III 

PCI 

PC II 

PC III 

GLS 

0.19 

0.45 

0.33 

-0.03 

0.57 

-0.11 

0.43 

-0.33 

0.15 

LNA 

-0.12 

-0.07 

0.57 

0.02 

0.30 

-0.24 

0.22 

-0.18 

-0.15 

LMX 

-0.07 

-0.19 

0.50 

0.01 

0.25 

-0.37 

0.22 

-0.15 

-0.01 

LFO 

-0.04 

0.44 

-0.22 

0.00 

0.26 

-0.01 

0.17 

-0.20 

0.53 

LZA 

0.00 

-0.15 

0.05 

0.03 

0.21 

0.16 

0.26 

0.02 

-0.05 
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Means, Males 

% Total: 

Age group 2 = 

74.2% 

Age group 3 = 89.4% 

Age group 4 = 

85.7% 

Var 

PC I 

PC II 

PC III 

PCI 

PC 11 

PC III 

PC I 

PC II 

PC III 

LBC 

0.20 

0.17 

-0.15 

-0.00 

0.42 

0.12 

0.22 

-0.12 

0.08 

WRO 

-0.03 

0.06 

0.16 

-0.02 

0.14 

-0.09 

0.12 

-0.07 

-0.20 

ICO 

0.02 

0.13 

0.03 

0.00 

0.03 

-0.03 

0.00 

0.00 

0.07 

A ZB 

0.39 

-0.15 

-0.14 

-0.03 

0.13 

0.02 

0.20 

-0.08 

-0.51 

SBP 

0.06 

0.21 

0.16 

0.01 

0.10 

-0.08 

0.07 

-0.02 

0.02 

LDI 

0.22 

-0.08 

0.18 

0.01 

0.14 

0.03 

0.17 

-0.08 

-0.30 

LMT 

0.07 

0.23 

0.04 

-0.00 

0.11 

0.04 

0.14 

-0.04 

0.18 

LPF 

-0.04 

0.12 

0.24 

-0.01 

0.11 

0.08 

0.13 

-0.02 

0.04 

BAM 

0.12 

0.33 

0.02 

-0.02 

0.02 

0.00 

0.01 

0.00 

0.41 

BAB 

0.36 

0.29 

0.05 

-0.01 

0.07 

0.08 

0.16 

-0.06 

-0.07 

HRO 

0.05 

0.11 

0.12 

-0.00 

0.08 

-0.03 

0.08 

-0.04 

0.06 

HMC 

0.15 

-0.13 

0.05 

-0.00 

0.18 

0.17 

0.19 

0.04 

-0.12 

MHS 

0.20 

0.03 

-0.06 

0.01 

0.09 

0.21 

0.13 

0.03 

-0.04 

HBP 

0.16 

0.05 

0.09 

-0.01 

-0.02 

0.15 

0.07 

-0.01 

0.16 

OLD 

0.42 

-0.21 

-0.06 

-0.12 

0.07 

0.74 

0.36 

-0.02 

-0.02 

LCP 

0.42 

-0.15 

0.06 

0.03 

0.25 

0.24 

0.36 

0.19 

0.05 

DHI 

0.05 

-0.03 

0.06 

0.23 

0.12 

-0.07 

0.08 

-0.00 

0.08 

AHM 

0.11 

0.03 

0.13 

0.52 

0.03 

-0.08 

0.15 

0.36 

0.09 

PHM 

0.24 

-0.22 

0.07 

0.81 

-0.03 

0.09 

0.27 

0.76 

0.03 

% PC 

45.1 

18.0 

11.1 

74.1 

9.4 

5.9 

69.7 

11.1 

4.9 


geographicaly closer GLs in some age and 
sex groups (M-X-A, A-X, M-I, M-I-D, O-E, 
V-M, P-N-L, P-N, P-L). From these, the most 
consistent were GLs A-X in all specimens, 
except for old males (age group 4). However, 
shorter phenetic distances were most com¬ 
mon among GLs geographically separated. 
Such groupings were not consistent in all age 
and sex groups (fig. 4). 

Cluster Analysis 

In the cluster analysis using only subadult 
(2) and adult (3) specimens from the same 
12 GLs (A, B, C, E, I, L, M, O, P, R, X, and 
Y), the localities clustered in several differ¬ 
ent subgroups according to age and sex. Such 
clusters were not related to the geographic 
nearness of the samples in every case (fig. 
5). The phenogram with the highest cophe- 
netic coefficient was that of subadult males 
(2, r = 0.915) followed by those of the fe¬ 
males of that age (2, r = 0.8728), adult fe¬ 
males (3, r = 0.8345), and adult males (3, r 
= 0.8062). 

Only the phenogram for subadult females 


could be compared to that of adult males (r 
=0.8519. fig. 5), since Mantel’s tests resulted 
in r < 0.8 in the other cases (Rohlf and Fish¬ 
er, 1968). In spite of this, certain common 
groupings were recognized among the four 
phenograms. The first case was a cluster in¬ 
cluding T. u. albigularis (A) and the pocket 
gophers from Sierra de Tlaxco (X) in adult 
females, as well as in the two male groups. 
Another group was formed by T. u. tolucae 
(O) with T. u. orizabae (R) in both adult fe¬ 
males and subadult males. Note that in every 
phenogram, T. u. pullus (P) separated either 
from the rest (all of the females and adult 
males) or from the annexed subgroup (sub¬ 
adult males). 

The clusters in figure 5 do not coincide 
completely with each other nor with the clus¬ 
ters in figure 4 depicted by discontinuous 
lines. In the latter, cluster analyses were run 
using all the GLs available for each age and 
sex group. Cophenetic coefficients (r) were 
0.93 for young females, 0.83 for adult fe¬ 
males, and 0.96 for old females. In males, r 
was 0.83 in young males, 0.90 in adult 
males, and 0.75 in old males. 
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THOMOMYS UMBRINUS 



Z = 44.05 

Z = 50.3 r = 0.17 

r = 0.68 p = 0.19 

p = 0.99 



Fig. 4. Distribution of populations of Thomomys umbrinus along the Mexican Transvolcanic Belt in 
multivariate space, according to sex (2/c?) and age (2-4). Notice that Thomomys umbrinus pullus, when 
present, always separated from other samples and that the position of the samples is variable in each 
sex and age. Axes and spheres indicate the position of group localities with respect to Principal Com¬ 
ponents Analyses computed on the population means. Lines connecting group localities result from a 
Minimum Spanning Tree projected onto the first principal components (I-III). Percentages refer to the 
amount of variation explained by each principal component. Discontinuous lines indicate groups ac¬ 
cording to Cluster Analyses developed in the sample means (r = cophenetic coefficient for the pheno- 
grams; n = sample size). Mantel’s Z tests (r = correlation; p = significance) and arrows indicate 
comparability between phenograms. See names of group localities (letters) in table 1. Diagrams were 
drawn with the same scale and are rotated to different degrees for the sake of clarity. 


For consistent groups, this cluster analysis 
coincided with the MST and PCA, except for 
some GLs that formed different arrange¬ 
ments. In young females, GL V clustered 
within the group formed by GLs C-L, N, R, 
X, and M, instead of linking with GLs Y and 
T, which in turn were clustered within dif¬ 
ferent groups. In adult females, GL T 
grouped with GLs I-M, O-R, L, and V, in¬ 
stead of grouping within the cluster including 
GLs A-X, C, S, B, and D. In old females, 
GLs 1 and V clustered closely together in¬ 
stead of gathering in separate groups (I, M, 
D, and V, O, R). In young males, the phen- 
ogram coincided perfectly with the MST. In 
adult males, GLs D-Y cluster closely instead 
of gathering with GLs A and E, respectively. 


In old males, the analyses were coincident, 
except for the shorter distance between GL 
D and GLs 1-M implied in the MST. 

Classification Discriminatory Analysis 

Table 8 shows the total number of speci¬ 
mens examined, as well as the percentages 
of correct classification according to their 
taxonomic designation or according the GL 
of origin. A total of 273 (74.2%, n = 368) 
males and 416 (71.1%, n = 585) females 
were correctly classified. All males from GLs 
P, L, T, D, Y, and S were completely classi¬ 
fied, as well as females from GLs, P, Y, and 
S. Percentage of incorrectly classified fe¬ 
males ranged from 50% in GL T to 4.4% in 
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p = 0.9279 

Fig. 5. Clustering of 12 populations of Thomomys umbrinus along the Mexican Transvolcanic Belt, 
according to sex (9 9/<3<3) and age (2, 3). Arrows show the comparisons of the phenograms using 
Mantel’s Z tests (r = correlation; p = significance). Cophenetic coefficients are r = 0.8728 in 92; r = 
0.8345 in 93; r = 0.9150 in <32; r = 0.8062 in <33. 


GL V, while males ranged from 46.5% in GL 
A to 12.4% in GL V. 

Most of the GLs with larger sample sizes 
had more degree of misclassification (X, A, 
I, and O in females; A, O, and X in males) 
and females had a greater number of wrongly 
classified individuals than males. Misclassi¬ 
fication in females occurred in 82.3% of all 
GLs (n = 14) and varied from one (GLs T 
and V) to 13 (GL X) with an average of 5.1 
GLs. GLs C, 1, O, A, and X had wrongly 
classified females in more than one-third of 
the total number of GLs. In males, misclas¬ 
sification occurred in 64.7% (n = 11) of the 
GLs, fluctuated from two (L and M) to 10 
(GL X) GLs, and averaged 5.2 GLs. GLs E, 
O, A, and X also had wrongly classified in¬ 
dividuals in more than one-third of the ex¬ 
amined populations. 

Similarly, more than being related to the 
geographic distribution of the GLs, the de¬ 


gree of overlapping of a certain GL is also 
related to sample size: the larger it is, the 
greater the overlapping (table 8). In general, 
except for the GLs P and D in males, GLs 
overlapped through poorly classified speci¬ 
mens. In the extremes, GL V (males, 3.4%; 
females, 4.4%) showed the lowest overlap¬ 
ping percentages and GL T (males, 90.0%; 
females, 66.7%) the highest. In females, 
overlapping of a GL with the others varied 
from one (GL P) to eight (GL O) but the 
average number was nine GLs (52.9% of the 
total number of GLs). Overlapping in males 
ranged from one (GLs V and M) to six (GL 
L) and the average was four GLs (23.5% of 
all GLs). 

Color 

Since color was not included in the anal¬ 
yses but has been cited frequently in the lit- 
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erature, some detailed notes on color are in¬ 
cluded here. In general, color does not cor¬ 
relate well with the taxonomic designations 
of specimens at UAMI. However, from a su¬ 
perficial qualitative analysis of the soil ex¬ 
tracted from the claws of some specimens, it 
was noted the hair shade might be related to 
the color of habitat soil. In areas where soil 
is yellowish since it contains more clay, 
specimens tend to be Clay 26, Tawny 38, 
Cinnamon 39, or Cinnamon Rufous 40. In 
darker soils with a larger amount of humus, 
specimens tend to be Vandyke Brown 121, 
221. In specimens from darker soils with vol¬ 
canic sand, specimens are also Vandyke 
Brown 121, and they can be Blackish Neutral 
Gray 82. All of the shades described are 
bright in specimens with shiny and, occa¬ 
sionally, greasy pelage but can be dull in oth¬ 
ers, especially in clay and sandy soil. 

Thomomys umbrinus albigularis (A). The 
holotype and topotypes show regular ochra- 
ceous coloration (Colors 26, 38, 39 and 40), 
on the whole back in the older specimens, 
but in the younger ones, the color shade 
tends to become darker toward the middle 
line (Color 121). Five of the six examined 
specimens showed a white tail tip. The backs 
of the hands and feet were white. The ventral 
coloration is lighter and blends with gray 
shades (like “dirty sand”). The hair base is 
Medium Plumbeus 87 and the distal half 
shows the described shade or it is white. The 
white spot on the throat and the internal skin 
of the cheek pouches are also white. The hair 
is silky and bright. 

Thomomys umbrinus martinensis (M). The 
cranium of the holotype had one printed label 
with the name of Thomomys umbrinus shel- 
doni and a handwritten label, by E. W. Nel¬ 
son, with the name of Thomomys umbrinus 
peregrinus. The skin is very deteriorated and 
the colors are mixed. In general, it is like that 
of T. u. albigularis, with molting indication 
in different brown shades. Coloration is ligh¬ 
ter (Cinnamon 39) on the coccyx area and on 
the flanks and darker (Tawny 38) on the 
head, back, and lumbar region. Because of 
the exposure of the hair base which is plum¬ 
beous gray (Color 87), the flanks show a 
light and grayish brown. From the waist to 
the anterior extreme and on the middle line, 
the holotype lacks hair. The tail and the ven¬ 


tral skin are as in T. a. albigularis and it 
shows a white stripe on the neck. The skin 
in the cheek pouches is whitish. It could 
seem either a young specimen from this sub¬ 
species which is molting or one of the light¬ 
est colored from T. u. peregrinus, T. u. to- 
lucae or 7. u. umbrinus, except that the color 
of the sides is brighter and more obvious. 
The other topotypical specimen is even 
younger and with molting evidence as well; 
thus coloration from its head down to the 
waist is dark and grayish. It has a very con¬ 
spicuous round, ocher brown (Color 39) spot 
on the lower half of the left leg. In general, 
the hair transition goes toward the reddish 
ochre (Color 38), lighter on the flanks. The 
ventral hair is like that of T. u. albigularis 
and white appears on the throat, inside the 
cheek pouches, and on the back of hands and 
feet. The hair is greasy. 

Thomomys umbrinus orizabae (R). The 
holotype is a melanistic specimen with a 
black matte shade (Color 82) and very dark 
brownish glittering areas on the head and 
middle line of the body (Color 121). It shows 
molting from the coccyx to the base of the 
tail where it has wooly juvenile hair. The 
dark color becomes lighter toward the flanks 
where it becomes ash black (Color 82). In 
contrast, half of the back of hands and feet 
as well as two-fifths of the tail are white. As 
in the flanks, the ventral pelage is also ash 
black (Color 82), but dark brown (Color 121) 
on the mammal gland areas and where the 
thighs begin. The mouth perimeter, the throat 
and the cheek pouches are white. Melanism 
is also present among the 14 topotypes, but 
two specimens, which show evidence of 
molting from head to waist, are ochre brown 
color (Color 26, 39), like T. u. albigularis or 
T. u. tolucae. All of the melanic specimens 
show a lighter hair base (Color 87) and dark¬ 
er at the tips (Color 82), as well as the con¬ 
trasting white color on the back of hands, 
feet and the tip of the tail. Five specimens 
show a small white spot on the throat, and, 
except for one, all have white on the inside 
skin of their cheek pouches. The ventral skin 
of the light specimens is as in T. u. albigu¬ 
laris. In the melanistic ones, it is plain gray 
(Dark Neutral Gray 83), and a little bit ligh¬ 
ter in both scapular and pelvic areas (Medi¬ 
um Neutral Gray 85). 
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TABLE 8 

Morphometric Overlap Among 17 Thomomys umbrinus Populations in the 
Mexican Transvolcanic Belt, According to Taxonomic Designation 

Group localities are mentioned from west to east (see table 1 for names). In the last two rows 1 indicates the total 
number of individuals considered as belonging to a certain taxon and tn refers to the total number of overlapping 
individuals for each case. The last column shows the number of specimens examined. Numbers within each cell 
indicate number of specimens and the percentage represented. Number of individuals and percentage of correct clas¬ 
sification in the original sample appear along the diagonal in bold face. 
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Thomomys umbrinus peregrinus (E). The 
holotype shows a dark brownish shade (Col¬ 
or 121), darker on the cheeks, nose, forehead 
and dorsum down to the waist, but especially 
on the back. It has a slight white spot above 
the right ear. The color gets lighter toward 
the flanks, cheekbones and below the waist. 
The tip of the tail, back of the hands and feet 


are white. It has a darker shade than T. u. 
albigularis and its belly is more grayish and 
whitish, though yellowish (Straw Yellow 56). 
Similar to that of T. u. albigularis, the brown 
back in the topotypes goes from ochre (Color 
38) to a dark shade (Color 121), but it is 
darker on the middle line. Moreover, the ab¬ 
domen is not quite different and 11 out of 
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TABLE 8 
(Continued) 
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the 15 specimens examined showed the 
white spot on the throat. The tip of the tail 
tends to be white, as well as the back of the 
hands and feet. The shade of this gopher is 
among the darker ones (Color 38) of T. u. 
albigularis, and the lightest ones (Color 26, 
39) of T. u. orizabae. Their hair is bright. 

Thomomys umbrinus pullus (P). In the 
case of the specimens from Los Tanques, 
similar patterns to those described above 
were found, but intermingled. Thus, some 
melanistic specimens (Color 82) were ob¬ 
served, and others were dark brown (Color 


121) as it was described for any of the sub¬ 
species with that shade. Their hair is bright. 

Thomomys umbrinus tolucae (O). The ho- 
lotype and topotypes show a general shade 
very similar to that of T. u. peregrinus, but 
with a duller shade (ash). The hair is greasy. 
The holotype shows a more grayish shade 
because it has lost a lot of hair and the gray 
(Color 84) base of hair is exposed at the ab¬ 
domen. In other ways, it would look like T. 
u. albigularis, to which it is also similar in 
the shade of the cheekbones, outside of the 
forearms and especially, on both sides of the 
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base of the tail. Five of the nine specimens 
examined show a white spot on their throat. 
In a few, the tail is white from the half and, 
in others, it is in the last third. Likewise, the 
skin on half the back of hands and feet and 
inside the cheek pouches is white. 

Thomomys umbrinus umbrinus (S). The 
color of the topotypes is quite variable, since 
most are young. The darkest ones are similar 
to T. u. tolucae with a dull-ash shade, and 
the lightest ones to T. u. albigularis. In all 
of them, the distal half of the tail is white 
and they show the white spot on their throat. 
The white back of hands and feet is more 
evident than in T. u. albigularis. The ventral 
shade is as in T. u. albigularis or as in T. u. 
peregrinus or T. a. tolucae. 

Thomomys umbrinus vulcanius (V). The 
holotype is dark ash brown (Color 121) due 
to the abundant gray undercoat hair (Color 
85) which gives it a wooly texture. The back 
is dark brown (Color 121) becoming lighter 
on the flanks. The tip of the tail and the back 
of hands and feet are white. The ventral color 
is as in Thomomys umbrinus peregrinus. The 
mouth perimeter and the skin of the chin are 
white. The hair is greasy. 

Thomomys umbrinus ssp. (C, X, I, T, B, 
Y). The same diversity of color patterns was 
found in the pocket gophers from elsewhere 
in the MTB and without taxonomic desig¬ 
nation. Melanistic (Color 82) specimens are 
also found in Thomomys umbrinus ssp. from 
Sierra de Las Cruces (C) and those from the 
Sierra de Tlaxco (X). Among the latter, there 
are also individuals with the ochre brownish 
types (Colors 26, 38, 39, 40) and dark brown 
(Color 121), which can be shiny or ash. Oth¬ 
er pocket gophers that show a darker brown 
(Color 121), as in T. u. orizabae, are those 
from Sierra de Tlaxco (X), Paso de Cortes 
(T. u. vulcanius, V), La Malinche (I), Totol- 
tepec (T), and those from the Cerro de San 
Andres (N) and Sierra de Las Cruces (C). 
Among those that have a brown ash shade 
are specimens from La Malinche (I), Aman- 
alco de Becerra (B), and Cerro de San An¬ 
dres. Individuals that show an ocher brown 
shade (Colors 26, 38, 39 and 40), such as T. 
u. albigularis, are from La Malinche (I), 
Oyameles (Y) and Sierra de Tlaxco (X). 


DISCUSSION 
Nongeographic Variation 

Genetic variation within a population can 
be seen as its fine tuning to habitat conditions 
in response to natural selection (Mayr, 1970). 
As a result, we assume that populations of 
the same species located in ecologically dis¬ 
tinct habitats will show nongeographical dif¬ 
ferences in ontogenetic variation, sexual di¬ 
morphism, and other individual variables 
such as coloration, shape, and size. There¬ 
fore, the analysis of nongeographic variation 
is a preamble to the study of geographic var¬ 
iation (Rarmrez-Pulido et al., 1991); but 
more importantly, it is one of the most ac¬ 
curate tools to understand the morphometric 
mechanisms involved in several aspects of 
the biology of populations (Owen and 
McBee, 1990; Patterson and Patton, 1990; 
Alves de Oliveira, 1992; Zelditch et al., 
1992; Lovich and Whitfield Gibbons, 1992). 

This is especially important in the case of 
Thomomys umbrinus (results shown here), T. 
bottae (Patton and Smith, 1990), and T. 
townsendii (Rogers, 1991a), since the high 
structural elasticity of their morphological 
features affects the expression of both onto¬ 
genetic variation and sexual dimorphism. 
Therefore, in these species as well as in oth¬ 
ers with the same characteristics, a poor 
knowledge of their intrapopulation variation 
could bias the results of interpopulation var¬ 
iation, and thus, make them useless. 

In general, individual variation of Thom¬ 
omys umbrinus from Sierra de Tlaxco, Tlax- 
cala (table 2) is similar to that found for other 
geomyids (Baker and Genoways, 1975; Hol¬ 
lander, 1990) and other rodents (Baumgard¬ 
ner and Schmidly, 1981; Stangl et al., 1991). 
However, CVs in this pocket gopher are of 
smaller magnitude and, thus, suitable for tax¬ 
onomic comparisons (Long, 1969). Like¬ 
wise, in this as in other species of geomyids 
(Grinnell, 1931; Davis, 1937, 1938rWilkins, 
1985; Rogers, 1991a) and rodents (Yancey et 
al., 1993), females show greater morpholog¬ 
ical uniformity than males. Eisenberg (1981) 
has suggested that morphological character¬ 
istics enable differential strategies between 
sexes related to the energy they devote to 
other biological processes (e.g., territoriality 
in males, gestation and nursing in females). 
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Age turned out to be the most important 
factor in nongeographic variation of Thom¬ 
omys umbrinus from Sierra de Tlaxco (tables 
3, 5, GL X). This pattern has also been found 
in Mus musculus (Straney, 1978) and other 
populations of T. umbrinus along the Mexi¬ 
can Transvolcanic Belt (unpublished data). In 
contrast. Smith and Patton (1988), Patton and 
Smith (1990), and Rogers (1991a) found that 
sexual dimorphism (tables 4, 5) is more rel¬ 
evant among adult and old adults in other 
species of Thomomys. It must be noted, 
though, that these authors a priori excluded 
younger classes from their studies, while we 
found separation among individuals from 
ages 2 to 4. Also, these authors followed a 
different statistical approach in order to in¬ 
crease sample sizes and their criteria for as¬ 
signing individuals to age groups (Daly and 
Patton, 1986) are different from ours (Castro- 
Campillo et al., 1993). 

With our approach, we found that both age 
and sex play an outstanding role in the non¬ 
geographic variation of Thomomys umbrinus 
from Sierra de Tlaxco. Indeed, they interact 
in one-third of the measurements (table 5). 
Therefore, it was decided to run the analyses 
on the geographic variation with separate age 
groups and sexes. Hollander (1990) arrived 
at a similar conclusion upon examining the 
variation of Cratogeomys castanops. 

Interaction between sex and age (table 5) 
suggests a differential pattern of develop¬ 
ment between males and females (Castro- 
Campillo et al., 1993). In Thomomys umbri¬ 
nus as in T. bottae (Patton and Smith, 1990), 
females seem to reach the proportions of the 
adult cranium sooner than males, which im¬ 
plies lower morphological variation among 
females than among males (Patton and 
Smith, 1990). This fact, which has been con¬ 
tinuously observed in other pocket gophers, 
has favored the use of females in compari¬ 
sons and descriptions of the taxa (Merriam, 
1895; Bailey, 1915; Grinnell, 1931; Davis, 
1937). 

The difference between skull growth of 
males and females of the congeneric species 
Thomomys bottae has been given various in¬ 
terpretation by different authors. To Howard 
and Childs (1959), it implied that males con¬ 
tinue to grow throughout their life. In con¬ 
trast, for Patton and co-workers (Daly and 


Patton, 1986; Patton and Smith, 1990), it 
proves the proportional increasing loss of 
smaller sized individuals within the popula¬ 
tion. In the case of Thomomys umbrinus, it 
is necessary to carry out growth studies to 
arrive at any conclusion. It will not be sur¬ 
prising if males continue to deposit calcium 
in their skulls at an age when females have 
finished doing so and are losing that element 
through pregnancy and lactation as in other 
mammals (Bentley, 1998). Moreover, sexual 
differences in skull conformation might be 
caused by less opportunity for females to de¬ 
posit calcium in their bones because of men¬ 
struation, as in human females. 

Finally, we must note that the skull devel¬ 
opment of Thomomys has been related to ex¬ 
trinsic patterns (Bailey, 1915; Patton and 
Brylski, 1987); indeed, the great elasticity 
shown by these gophers is multifactorial 
(Smith and Patton, 1988; Lessa and Patton, 
1989). Therefore, both ontogenetic variation 
and secondary sexual differentiation should 
be examined in depth, as done here, before 
studying geographic variation. Such scope 
will allow immediate factors, such as the en¬ 
vironmental influence, to be considered sep¬ 
arately from nonimmediate factors, such as 
evolutionary and adaptive processes devel¬ 
oped throughout the phylogenetic history of 
these species (Patton and Smith, 1990). 

Geographic Variation 

Working with subspecies, and especially 
using morphometry, involves a careful inter¬ 
pretation of the existence of gene flow on the 
basis of phenetic similarity (Rogers, 1991a). 
Specimens assigned to each subspecies 
should be reasonably uniform (Pimentel, 
1979; Manly, 1986). This supposition is not 
fulfilled in the analyzed populations, partly 
because of the different sample sizes. 

The diversity of results in the multivariate 
analyses (figs. 3—5, tables 6, 7) may be 
caused, in part, by differing sample sizes, 
very small in some GLs (Oyameles, Totol- 
tepec, T. u. umbrinus). However, it was de¬ 
cided to use all of the specimens to explore 
the phenetic similarity among the popula¬ 
tions. Other authors (Smith et al., 1983; 
Smith and Patton, 1988; Hollander, 1990) 
have also included single specimens or very 
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small sample sizes (n < 4) in analyses of 
patterns among geomyids, transforming the 
data into their logarithmic expression (Neff 
and Marcus, 1980; Rohlf and Bookstein, 
1990). Homogenization of quantities is a use¬ 
ful strategy in cases where “all of the adult 
specimens” have been used. However, since 
adult specimens were separated here into six 
groups, according to age and sex classes, var¬ 
iability is divided and the results are less ex¬ 
treme (R. D. Owen, and R. O. Gonzalez-Ro- 
bles, personal commun.). 

Likewise, the use of six age groups re¬ 
sulted in a better expression of their geo¬ 
graphic variation (i.e., high cophenetic co¬ 
efficients in the cluster analysis or in the ei¬ 
genvalues of the PC A, table 7). Moreover, 
splitting of variation indicates that in the 
PCA on the GL means (fig. 4, table 7), PC I 
was a factor of shape (i.e., eigenvalues of 
different signs and some variables distinctly 
heavier than others; Pimentel, 1979; Neff and 
Marcus, 1980) and not exclusively of size, as 
is usually the case. 

Results of classification discriminant anal¬ 
yses of variables (table 7) confirm the gen¬ 
eral pattern found in the rest of the analyses 
(tables 6—7, figs. 3—5) in terms of a general 
overlap among populations. Indeed, although 
discrimination accuracy was good (> 70%), 
population sample size had a more definite 
influence, since all of the GLs with n < 10, 
except for Totoltepec (GL T), are better clas¬ 
sified. However, these specimens, except for 
T. u. pullus, also overlapped more with other 
larger samples. Indeed, the greater the sam¬ 
ple size, the larger the overlap with other 
small and large samples. 

The fact that in Thomomys umbrinus, male 
dimensions from different localities overlap 
less than female dimensions, confirms that 
the latter are morphometricaly more uniform 
(Castro-Campillo et al., 1993) among GLs 
(table 8). This is a phenomenon also dem¬ 
onstrated in other species of the genus (Grin- 
nell, 1931; Davis, 1937; Patton and Smith, 
1990; Rogers, 1991a) and in other genera of 
pocket gophers (Davis, 1940; Honeycutt and 
Schmidly, 1979; Wilkins, 1985; Hollander, 
1990). As mentioned before, this can be re¬ 
lated to different sexual roles (Eisenberg, 
1981; Castro-Campillo, et al., 1993). To con¬ 
firm this, it will be necessary to carry out 


both ethological and ecological studies in 
situ, as well as physiological studies. 

Overall, with the exception of T. u. pullus, 
results of multivariate analyses also show no 
definite pattern attributable to geographic 
distribution (figs. 3—, table 6). Different GLs 
are grouped together regardless of geograph¬ 
ic distance between them according to the 
sex and age considered in the analysis. The 
only taxon consistently segregated in every 
analysis, when it is present, is T. u. pullus. 

The absence of a common pattern and the 
diversity of combinations of GLs in multi¬ 
variate space (figs. 3—5) may suggest differ¬ 
ent factors in the association of individuals 
in the different ontogenetic and sexual 
groups. It also may suggest that chance and 
small sample sizes are important. These re¬ 
sults can be interpreted as a structural ho¬ 
mogeneity referable to shape and size. Size 
is the most important characteristic to distin¬ 
guishing T. u. pullus from the rest. 

Rogers (1991a) also found morphometric 
homogeneity in Thomomys townsendii, in 
which she examined nine subspecies in four 
states of the United States. Since size was 
very important in her analysis, she decided 
to interpret her results with caution because 
Patton and co-workers have shown that in 
Thomomys bottae, size depends on ecology 
(Straney and Patton, 1980; Smith et al., 
1983; Kennedy and Lindsay, 1984; Daly and 
Patton, 1986; Patton and Brylski, 1987; Pat¬ 
ton and Smith, 1989, 1990). After examining 
the genic variation of Thomomys townsendii, 
Rogers (1991b) concluded that there are two 
and not nine subspecies. 

After examining 43 subspecies of Thom¬ 
omys bottae in California, Patton and Smith 
(1990) found that cranial size is a character 
that responds in a very plastic way, and that 
therefore, shows great variation. They inter¬ 
preted this to mean that size is related to non- 
genetic effects such as nutrition and other en¬ 
vironmental factors. In contrast, shape ex¬ 
presses allometric differences in the mensu¬ 
rable characters; according to Smith and 
Patton (1988) it is more useful in defining 
regional geographic units, even though cer¬ 
tain specific morphologies are recurrent 
along the species range. From biochemical 
analyses, these authors have found an inverse 
relationship between the magnitude of ge- 
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netic distance and geographic distance 
(Smith et al., 1983). Based on their findings, 
and with a conceptual modification of Grin- 
nell’s criteria (1935) for the recognition of 
the subspecies. Smith et al. (1983) recog¬ 
nized 15 infraspecific entities of Thomomys 
bottae in California. 

In the case of Thomomys umbrinus, the ex¬ 
tensive overlap revealed in the ANOVA and 
most of the multivariate analyses (figs. 3, 4, 
table 8), demonstrates that skull shape of this 
species exhibits no clear geographic pattern 
along the MTB. The only exception is T. u. 
pulliis. 

Davis (1938) noted that the skull of Thom¬ 
omys bottae and Thomomys quadrattus (a cur¬ 
rent subspecies of T. talpoides ) is larger in 
populations inhabiting lower altitudes and 
deeper soils. In contrast, the skull is smaller 
in those populations found in high altitudes, 
where the ground is shallower (Davis, 1938). 
Besides size, Davis (1938) mentioned the 
same relationship in the presence and con¬ 
spicuousness of the skull crests (sagittal, 
lambdoid, and temporal crests), and he con¬ 
sidered it inappropriate to designate a trino¬ 
mial on the basis of these crests. Davis (1938) 
concluded that variation within the species ex¬ 
amined, together with the recognition of the 
topography and geological history, was suf¬ 
ficient to determine the validity of a subspe¬ 
cific entity. Choate and Williams (1978) 
reached a similar conclusion in the interpre¬ 
tation of variation in Microtus ochrogaster. 

Davis’s conclusions (1938) make sense for 
Thomomys umbrinus pullus, since the skull is 
smoother and smaller in this subspecies, a fact 
possibly related to the depth of the soil and 
to a less intricate topography of its habitat 
(Pasquare et al., 1987a, 1987b). To a lesser 
degree, the differences found in T. u. pullus 
can also be related to the geological history 
of its distribution area (Demant, 1978, 1982; 
Nixon et al., 1987; Pasquare et al., 1987b). 

Further evidence of cranial homogeneity in 
Thomomys umbrinus from the Mexican Trans- 
volcanic Belt is that all populations located 
east of Patzcuaro, Michoacan, show recurrent 
shapes even when they are separated by more 
than 100 km. Such is the case of populations 
from Amanalco de Becerra (GL B) and La 
Malinche (GL I). This has been interpreted as 
canalization within a fossorial niche (Honey¬ 


cutt and Schmidly, 1979; Lessa and Patton, 

1989) and, thus, is given less weight than the 
allometric patterns in defining infraspecific 
entities (Patton and Smith, 1990). 

With respect to pelage coloration, there is 
overlap regardless of taxonomic designation 
or geographic location. Moreover, examina¬ 
tion of the holotypes and topotypes makes 
clear that color discrepancies, rather than be¬ 
ing diagnostic features, are conditions related 
to age of the specimen, the presence of molt¬ 
ing, and the small number of specimens ex¬ 
amined. Interestingly, the larger the sample 
size, the more color types are found: as in 
pocket gophers from Sierra de Tlaxco (GL 
X), T. u. albigularis (A), and T. u. pullus (P). 
Color patterns were not considered in the 
original descriptions or in the group reviews 
(Bailey, 1915; Nelson and Goldman, 1934; 
Hall and Villa, 1948). In Thomomys umbri¬ 
nus, as in other pocket gopher species, differ¬ 
ent color types are not useful in characterizing 
taxonomic entities, which may share different 
color types in the MTB. Therefore, coloration 
in such species may not be a useful diagnostic 
feature; rather, it is related to the color and 
content of the soil (Davis, 1937; Rogers, 
1991a; Patton and Smith, 1990). 

In summary, results support the distinction 
of two subspecies (sensu Grinnell, 1935; 
Mayr, 1969) of Thomomys umbrinus on the 
southern boundary of its distribution. One of 
these, Thomomys umbrinus pullus, includes 
the populations in the surroundings of Patz¬ 
cuaro, Michoacan. The other includes all of 
the other populations eastward from Sierra 
de Las Cruces, Michoacan, to the Pico de 
Orizaba Volcano, Veracruz. This second tax¬ 
on, which includes several subspecies for¬ 
merly recognized, as well as populations 
without previous designation, should be re¬ 
ferred to Thomomys umbrinus umbrinus, the 
oldest name available. However, it is neces¬ 
sary to check the situation of this name, since 
the designation of the specimens of Boca del 
Monte is not very clear (Bailey, 1906), and 
the neotype(?) assigned to it was lost (D. 
Wilson, personal commun., November 

1990) . 

Acknowledgment of two subspecies in¬ 
stead of eight in the Mexican Transvolcanic 
Belt agrees in general with the biochemical 
and cytological evidence found by Hafner et 
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al. (1987). They showed that the Thomomys 
umbrinus specimens examined in the Mexi¬ 
can Transvolcanic Belt belong to the same 
cy to type (FN = 78) and that the genic dif¬ 
ferences between them are minimal (in fact, 
they place them in the same electrophoretic 
group). Likewise, Hellenthal and Price 
(1984) have found specific associations of 
ectoparasites shared by all of the Thomomys 
umbrinus specimens examined in the Mexi¬ 
can Transvolcanic Belt, as opposed to those 
located a few kilometers north. 

Perhaps relationships found among popu¬ 
lations of Thomomys umbrinus in the Mexi¬ 
can Transvolcanic Belt can be explained by 
the geological history of the Quaternary Age 
of the area, which also explains populations 
of other species of the genus (Davis, 1937; 
Rogers, 1991b). While the formation of the 
Mexican Transvolcanic Belt goes beyond the 
late Miocene (Cebull and Shurbet, 1987; 
Nixon et al, 1987), the origin of vulcanism 
in the region began during the Pliocene (De- 
mant, 1978; 1982; Nixon et al., 1987). On 
the other hand, it is known that the Geomyi- 
dae originated in the late Miocene (Kurten 
and Anderson, 1980), but the oldest records 
of the species in Mexico are referable to the 
Late Pleistocene (Barrios Rivera, 1985). 
Some of these fossil records are found within 
the study area at Tequixquiac (Hibbard, 
1955) and Tequesquinahua (Alvarez, 1966), 
both in the State of Mexico. 

In retrospect, one can visualize a large 
zone along the valleys of Mexico, Toluca, 
and Pachuca having a homogeneous environ¬ 
ment, with both climate and vegetation sim¬ 
ilar to those inhabited by the recent pocket 
gophers. This area was then reduced and par¬ 
titioned by the changes in climate during the 
Quaternary Period (Hibbard, 1955). More¬ 
over, at the end of the period, the rise of the 
strato-volcanoes in the central and eastern re¬ 
gion of the Belt, less than one million years 
ago (Demant, 1978, 1982), separated Thom¬ 
omys umbrinus populations. 

Evidence of a larger duration of the cli¬ 
matic and ecological conditions is the fossil- 
iferous locality of Tequesquinahua, a site 3.2 
km northeast of Tlalnepantla (Alvarez, 
1966). Even though xeric species of plants 
are currently found there, in the past, Te¬ 
quesquinahua may have been covered by co¬ 
niferous forests and grassland, because rem¬ 


nants of Thomomys umbrinus and Neotomo- 
don alstoni were found there (Alvarez, 
1966). Also, the fossil site at Tequixquiac, 
near Zumpango in the State of Mexico (Hib¬ 
bard, 1955), has similar characteristics to Te¬ 
quesquinahua; specimens of Thomomys um¬ 
brinus from the Pleistocene have been re¬ 
ported there (Hibbard, 1955). 

Differentiation of Thomomys umbrinus 
pullus could be related either to this overall 
picture or to the nature of the soil that it in¬ 
habits. That is, this subspecies lives in a re¬ 
gion with less deposition of volcanic sand 
than in areas of other populations (Pasquare 
et al., 1987a, 1987b). 

CONCLUSIONS 

Both sexual dimorphism and, especially, 
age variation are important determining fac¬ 
tors in the variation of Thomomys umbrinus. 
Such variation seems to be linked to a very 
fine tuning between ontogenetic development 
and sexual roles. 

Analyses of geographic variation demon¬ 
strate that, except for Thomomys umbrinus 
pullus, previously recognized subspecies and 
other populations in the MTB are morpho¬ 
metrically similar regardless of their geo¬ 
graphic origin. This finding does not support 
systematic recognition of eight subspecies in 
that region, despite some degree of geo¬ 
graphic isolation. Besides, isolation among 
type localities is less than previously sup¬ 
posed, as shown by specimens obtained from 
intermediate localities. 

Hence, we recognize Thomomys umbrinus 
pullus and the nominal subspecies as the only 
two subspecies of Thomomys umbrinus oc¬ 
curring in the MTB. The nomenclature may 
be summarized as follows: 

Thomomys umbrinus umbrinus (Richardson) 

1829. Geomys umbrinus Richardson, Fauna 
Boreali-Americana 1: 202, type locality 
undefined. Bailey (1906: 4) restricted 
the type locality to southern Mexico, 
probably in the vicinity of Boca del 
Monte, Veracruz. General distribu¬ 
tion: Veracmz. 

1893. Thomomys umbrinus orizabae Merri- 
am, Proc. Biol. Soc. Washington 8: 
145, type locality: Mt. Orizaba, 9,500 
ft, Puebla. General distribution: 
Puebla. 
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1893. Thomomys umbrinus peregrinus Mer- 
riam, Proc. Biol. Soc. Washington 8: 
146, type locality: Salazar, 10,300 ft, 
Mexico. General distribution: Distri¬ 
to Federal, Hidalgo, Mexico, and Mo¬ 
relos. 

1934. Thomomys umbrinus albigularis Nel¬ 
son and Goldman, J. Mammal. 15: 106, 
type locality: El Chico, 9,800 ft. Sierra 
de Pachuca, Hidalgo. General distri¬ 
bution: Hidalgo, Puebla, and Veracruz. 

1934. Thomomys umbrinus martinensis Nel¬ 
son and Goldman, J. Mammal. 15: 108, 
type locality: San Martin Texmelucan, 
7,400 ft, Puebla. General distribu¬ 
tion: Known only from type locality. 

1934. Thomomys umbrinus tolucae Nelson 
and Goldman, J. Mammal. 15: 109, 
type locality: Northern slope of Toluca 
Volcano, 9,500 ft, Mexico. General 
distribution: Mexico. 

1934. Thomomys umbrinus vulcanius Nelson 
and Goldman, J. Mammal. 15: 109, 
type locality: Popocatepetl, 12,900 ft, 
Mexico. General distribution: Mexi¬ 
co, Morelos, and Puebla. 

Thomomys umbrinus pullus Hall and Villa 

1948. Thomomys umbrinus pullus Hall and 
Villa, Univ. Kansas Publ. Mus. Nat. 
Hist. 1: 251, type locality: 5 mi S Patz- 
cuaro, 7,800 ft, Michoacan. General 
distribution: Known only from the vi¬ 
cinity of type locality. 

ACKNOWLEDGMENTS 

Several colleagues and friends participated 
in some phase of this research. To all of 
them, we extend our most sincere and deep 
appreciation. 

Dr. Michael Smolen (Texas A&M Univer¬ 
sity) made valuable suggestions at the begin¬ 
ning of this study. Dr. Robert D. Owen gen¬ 
erously allowed ACC to compute the statis¬ 
tical analyses in his Texas Tech Univ. labo¬ 
ratory and made interesting comments to the 
research. M. S. Rosa Obdulia Gonzalez-Ro- 
bles also made helpful comments on the in¬ 
terpretation of the statistical results. 

Dr. Margaret Smith (UC Berkeley) helped 
to clear up some points and M. S. Ticul Al¬ 
varez (ENCB, IPN), brought out some as¬ 
pects that permitted a better understanding of 
the problem. Drs. Joaquin Arroyo-Cabrales 
(IN AH), Bernardo Villa (IB UN AM), Ismael 
Ferrusqufa-Villafranca (Inst. Geol., UNAM), 


Mary Ellen Holden (AMNH), and Sydney 
Anderson (AMNH) gave us some literature 
references that enriched the study. We thank 
Eric M. Brothers for editing the figures. A 
thorough revision of the manuscript by Drs. 
James Patton and Sydney Anderson together 
with a third anonymous reviewer undoubt¬ 
edly contributed to a more succinct and more 
understandable final version of the paper. 

A careful examination of specimens was 
made possible by the kindness and hospital¬ 
ity of several curators: M. S. Ticul Alvarez 
at CB; Drs. Don Wilson, Michael Carleton, 
Richard Thorington, and other members of 
the staff at USNM; Drs. Guy Musser and 
Sydney Anderson at AMNH; Dr. Robert 
Timm at KU; Dr. George Baumgardner at 
TCWC, and Dr. Stephen L. Williams at TTU. 

Technicians Juan Patino-Rodrfguez and 
Benjamin Vieyra-Rosas collected the speci¬ 
mens in the field and were helped by several 
students during the three years of fieldwork. 
Hugo Martinez-Paz, Elsa Gonzalez, Beatriz 
Gutierrez, and Claudia Aguilar helped us 
with word-processing of the Spanish version 
of the manuscript. Mr. Fernando Vega made 
the figures. Mr. Antonio H. Cabrera took care 
of an early draft of English version. 

This research was conducted to satisfy the 
Ph.D. dissertation requirements of ACC. It 
was financed by Consejo Nacional de Cien- 
cia y Tecnologfa (1253-N9203 and 400200- 
5-R29117N) and Secretarfa de Educacion 
Public a (93-01-09-002-177, 94-01-09-002- 
248). 

REFERENCES 

Alvarez, T. 

1966. Roedores fosiles del Pleistoceno de Te- 
quesquinahua, Estado de Mexico, Mex¬ 
ico. Acta Zool. Mex. 8: 1—16. 

Alves De Oliveira, J. 

1992. Estrutura da variagao craniana em po¬ 
pulates de Bolomys lasiurus (Lund, 
1841) (Rodentia: Cricetinae) do nord- 
este do Brasil. Tesis de Maestrla. Univ. 
Federal Do Rio de Janeiro, Rio de Ja¬ 
neiro, Brasil, 107 pp. 

Anderson, S. 

1966. Taxonomy of gophers, especially 
Thomomys in Chihuahua. Mexico. Syst. 
Zool. 15: 189-198. 

Bailey, V. 

1906. Identity of Thomomys umbrinus (Rich- 



34 


AMERICAN MUSEUM NOVITATES 


NO. 3297 


ardson). Proc. Biol. Soc. Washington 
19: 3-6. 

1915. Revision of the pocket gophers of the 
genus Thomomys. N. Am. Fauna 39: 1— 
136. 

Baker, R. J., and H. H. Genoways 

1975. A new subspecies of Geomys bursarius 
(Mammalia: Geomyidae) from Texas 
and New Mexico. Occas. Pap. Mus. 
Texas Tech Univ. 29: 1-18. 

Barrios Rivera, H. 

1985. Estudio analltico del registro paleover- 
tebradologico de Mexico. Tesis de Li- 
cenciatura no publicada. Univ. Nacion- 
al Autonoma de Mexico, Mexico, D. F., 
474 pp. 

Baumgardner, G. D.. and D. J. Schmidly 

1981. Systematics of the southern races of 
two species of kangaroo rats (Dipodo- 
mys compactus and D. ordii). Occas. 
Pap. Mus. Texas Tech Univ. 73: 1-27. 

Bentley, P. J. 

1998. Comparative vertebrate endocrinology. 
3 rd ed. Cambridge, Cambridge Univ. 
Press, 526 pp. 

Brown, W. L., and E. O. Wilson 

1954. The case against the trinomen. Syst. 
Zool. 3: 174-176. 

Castro-Campillo, A. 

1995. Taxonomfa de las pequenas tuzas 
Thomomys umbrinus en el Eje Volcan- 
ico Transverso. Ph. D. Thesis, Facultad 
de Ciencias, Univ. Nacional Autonoma 
de Mexico. Mexico, D. F.. 100 pp. 

Castro-Campillo, A., O. Gonzalez-Robles, and J. 
Ramlrez-Pulido 

1993. Determination of age groups in Thom¬ 
omys umbrinus (Rodentia: Geomyi¬ 
dae). Texas J. Sci. 45: 217-230. 

Cebull, S. E„ and D. H. Shurbet 

1987. Mexican Volcanic Belt: an interplate 
transform?. Geofs. Int. Spec. vol. on 
Mexican Volcanic Belt—Pt. 3A (S. P. 
Verma, ed.), vol. 26: 1-13. 

Choate, J. R., and S. L. Williams 

1978. Biogeographic interpretation of varia¬ 
tion within and among populations of 
the prairie vole, Microtus ochrogaster. 
Occas. Pap. Mus. Texas Tech Univ. 49: 
1-25. 

Daly, J. C., and J. L. Patton 

1986. Growth, reproduction, and sexual di¬ 
morphism in Thomomys bottae pocket 
gophers. J. Mammal. 67: 256-265. 

Davis, W. B. 

1937. Variations in Townsend pocket go¬ 
phers. J. Mammal. 18: 145—158. 

1938. Relation of size of pocket gophers to 
soil and altitude. Ibid. 19: 338-342. 


1940. Distribution and variation of pocket go¬ 
phers (genus Geomys ) in the south¬ 
western United States. Texas Agric. 
Exp. Stn. Bull. 590: 1-38. 

Demant, A. 

1978. Caracterfsticas del Eje Neovolcanico 
Transmexicano y sus problemas de in¬ 
terpretation. Univ. Nac. Auton. Mexico 
Inst. Geol. Rev. 2: 172-187. 

1982. Interpretation geodinamica del volcan- 
isrno en el Eje Neovolcanico Trans¬ 
mexicano. Ibid. 5: 217-222. 

Eisenberg, J. E. 

1981. The mammalian radiations. An analysis 
of trends in evolution, adaptation, and 
behavior. Chicago: Univ. Chicago 
Press, 610 pp. 

Fa, J., and L. M. Morales 

1991. Mammals and protected areas in the 
Trans-Mexican Neovolcanic Belt. In M. 
Mares and D. J. Schmidly (eds.). Top¬ 
ics in Latinamerican Mammalogy: his¬ 
tory, biodiversity, and conservation: 
199—226. Norman: Univ. Oklahoma 
Press, 468 pp. 

Garcia, E. 

1981. Modificaciones al sistema de clasifica- 
con climatica de Koppen (para adaptar- 
lo a las condiciones de la Republica 
Mexicana. Mexico, D. F.: Inst. Geogr., 
Univ. Nac. Auton. Mex., 252 pp. 

Gower, J. C., and G. J. S-Ross 

1969. Minimum spanning trees and single¬ 
linkage cluster analysis. Appl. Stat. 18: 
54-64. 

Grinnell, J. 

1931. A new pocket gopher from southeast¬ 
ern California. Univ. California Publ. 
Zool. 38: 1-10. 

1935. Differentiation in pocket gophers of the 
Thomomys bottae group in California 
and southern Oregon. Ibid. 40: 403- 
416. 

Hafner, M. S., J. C. Hafner. J. L. Patton, and M. 
F. Smith 

1987. Macrogeographic patterns of genetic 
differentiation in the pocket gopher 
Thomomys umbrinus. Syst. Zool. 36: 
18-34. 

Hall, E. R. 

1981. The mammals of North America, vol. 
1. New York: Wiley. 

Hall, E. R„ and K. R. Kelson 

1959. The mammals of North America. Vol. 
1. New York: Ronald Press. 

Hall, E. R„ and W. C. Russell 

1933. Dermestid beetles as an aid in cleaning 
bones. J. Mammal. 14: 372-374. 



2000 


CASTRO-C. AND RAMIREZ-P.: THOMOMYS UMBRINUS IN MTB 


35 


Hall, E. R„ and B. Villa R. 

1948. A new pocket gopher ( Thomomys ) and 
a new spiny pocket mouse ( Liomys) 
from Michoacan, Mexico. Univ. Kan¬ 
sas Publ. Mus. Nat. Hist. 1: 249-256. 

Heaney, L. R.. and R. M. Timm 

1983. Relationships of pocket gophers of the 
genus Geomys from central and north¬ 
ern Great Plains. Univ. Kansas Mus. 
Nat. Hist. Misc. Publ. 74: 1-59. 

Hellenthal, R. A., and R. D. Price 

1984. Distributional associations among Geo- 
mydoecus and Thomomydoecus lice 
(Mallophaga: Trichodectidae) and 
pocket gopher hosts of the Thomomys 
bottae group (Rodentia: Geomyidae). J. 
Med. Entomol. 21: 432-446. 

Hibbard, C. W. 

1955. Pleistocene vertebrates from the Upper 
Becerra (Becerra Superior) Formation, 
Valley of Tequisquisc, Mexico, with 
notes on other Pleistocen forms. Con- 
trib. Mus. Paleontol. Univ. Michigan 
12: 47-96. 

Hoffmeister, D. F. 

1969. The species problem in Thomomys bot- 
tae-Thomomys umbrinus complex of 
pocket gophers in Arizona. In J. K. 
Jones, Jr. (ed.). Contributions in mam¬ 
malogy: 75-91. Univ. Kansas Mus. 
Nat. Hist. Misc. Publ. 51: 428 pp. 

Hollander, R. R. 

1990. Biosystematics of the Yellow-eared 
pocket gopher, Cratogeomys castanops 
(Rodentia: Geomyidae) in the United 
States. Spec. Publ. Mus. Texas Tech 
Univ. 33: 1-62. 

Honeycutt, R. L., and D. J. Schmidly 

1979. Chromosomal and morphological vari¬ 
ation in the Plains Pocket Gopher. Geo¬ 
mys bursarius , in Texas and adjacent 
states. Occas. Pap. Mus. Texas Tech 
Univ. 58: 1-54. 

Howard, W. E„ and H. E. Childs. Jr. 

1959. Ecology of pocket gophers with em¬ 
phasis on Thomomys bottae mewa. Hil- 
gardia 29: 277-358. 

Kennedy, M. L.. and S. L. Lindsay 

1984. Morphologic variation in the raccoon, 
Procyon lotor , and its relationship to 
genic and environmental variation. J. 
Mammal. 65: 195—205. 

Kurten, B., and E. Anderson 

1980. Pleistocene mammals of North Ameri¬ 
ca. New York: Columbia Univ. Press, 
442 pp. 

Lapointe, F. J., and P. Legendre 

1992. Statistical significance of the matrix 


correlation coefficient for comparing 
independent phylogenetic trees. Syst. 
Biol. 41: 378-384." 

Leamy, L. 

1983. Variance partitioning and effects of sex 
and age on morphometric traits in ran- 
dombred house mice. J. Mammal. 64: 
55-61. 

Lessa, E., and J. L. Patton 

1989. Structural constrains, recurrent shapes, 
and allometry in pocket gophers (genus 
Thomomys). Biol. J. Linnean Soc. 36: 
349-363. 

Long, C. A. 

1969. An analysis of patterns of variation in 
some representative Mammalia. In J. 
K. Jones, Jr. (ed.). Contributions in 
mammalogy. Pt. 2. Studies on the na¬ 
ture and correlation of measures of var¬ 
iation: 89—302. Univ. Kansas Mus. 
Nat. Hist. Misc. Publ. 51: 428 pp. 

Lovich. J. E., and J. Whitfield Gibbons 

1992. A review of techniques for quantifying 
sexual size dimoiphism. Growth Dev. 
Aging 56: 269-281. 

Manly, B. F. J. 

1986. Multivariate statistical methods: a 
primer. New York: Chapman and Hall, 
159 pp. 

Mantel. N. A. 

1967. The detections of disease clustering and 
a generalized regression approach. 
Cancer Res. 27: 209—220. 

Mayr, E. 

1969. Principles of systematic zoology. New 
York: Mac Graw-Hill, 471 pp. 

1970. Populations, species, and evolution. An 
abridgment of animal species and evo¬ 
lution, 2 nd ed. Cambridge, MA: Belk¬ 
nap Press, 453 pp. 

Mayr, E.. and P. D. Ashlock 

1991. Principles of systematic zoology. 2 nd 
ed. New York: McGraw-Hill. 475 pp. 

Merriam, C. H. 

1893. Preliminary description of four new 
mammals from southern Mexico, col¬ 
lected by E. W. Nelson. Proc. Biol. Soc. 
Washington 8: 143—146. 

1895. Monographic revision of the pocket go¬ 
phers, Family Geomyidae (excl. of the 
species of Thomomys). N. Am. Fauna 
8: 258 pp. 

Neff, N. A., and L. F. Marcus 

1980. A survey of multivariate methods for 
Systematics. For a workshop. Numeri¬ 
cal methods in systematic mammalogy. 
Shippensburg, PA: Am. Soc. Mammal., 
243 pp. 



36 


AMERICAN MUSEUM NOVITATES 


NO. 3297 


Nelson, E. W., and E. A. Goldman 

1934. Pocket gophers of the genus Thomomys 
of Mexican mainland and bordering 
territory. J. Mammal. 15: 105-124. 

Nixon, G. T., A. Demant, R. L. Armstrong, and J. 
E. Harakal 

1987. K-Ar and geologic data bearing on the 
age and evolution of the Transmexican 
Volcanic Belt. In S. P. Verma (ed.), 
Geofs. Int. Spec. vol. on Mexican Vol¬ 
canic Belt—Pt. 3, vol. 26: 109-158. 

Ott, L. 

1984. An introduction to statistical methods. 
2 nd ed. Boston, MA: Duxbury Press, 
775 pp. 

Owen, R. D., and K. Mcbee 

1990. Analysis of asymmetry and morpho¬ 
metric variation in natural populations 
of chromosomedamaged mice. Texas J. 
Sci. 42: 319-322. 

Pasquare, G., L Vezzoli, and A. Zanchi 

1987a. Morphological and structural model of 
Mexican Volcanic Belt. In S. P. Verma 
(ed.), Geofs. Int. Spec. vol. on Mexican 
Volcanic Belt—Pt. 3B, vol. 26: 159— 
176. 

Pasquare, G., L Ferrari. V. Perazzoli, M. Tiberi, 
and F. Turchetti 

1987b. Morphological and structural analysis 
of the central sector of the Transmexi¬ 
can Volcanic Belt. In S. P. Verma (ed.), 
Geofs. Int. Spec. vol. on Mexican Vol¬ 
canic Belt—Pt. 3B, vol. 26: 177-193 

Patterson, B. D., and J. L. Patton 

1990. Fluctuating asymmetry and allozymic 
heterozygosity among natural popula¬ 
tions of pocket gophers {Thomomys 
bottae). Biol. J. Linnean Soc. 40: 21—36. 

Patton. J. L. 

1973. An analysis of natural hybridization be¬ 
tween the pocket gophers Thomomys 
bottae and Thomomys umbrinus in Ar¬ 
izona. J. Mammal. 54: 561—584. 

1993. Family Geomyidae. In D. E. Wilson 
and D. A. M. Reeder (eds.). Mammal 
species of the world. A taxonomic and 
geographic reference. 2 nd ed.: 469-476. 
Washington, D. C.: Smithsonian Insti¬ 
tution Press, xviii + 1206 pp. 

Patton, J. L., and P. V. Brylski 
1987. Pocket gophers in alfalfa fields: causes 
and consequences of habitat-related 
body size variation. Am. Nat. 130: 
493-506. 

Patton. J. L., and R. E. Dingman 

1968. Chromosome studies of pocket gophers 
genus Thomomys I. The specific status 
of Thomomys umbrinus (Richardson) in 
Arizona. J. Mammal. 49: 1-13. 


Patton. J. L., and J. H. Feder 

1978. Genetic divergence between popula¬ 
tions of the pocket gopher, Thomomys 
umbrinus (Richardson). Z. Saeuge- 
tierkd. 43: 17-30. 

Patton. J. L., and M. F. Smith 

1989. Population structure and genetic and 
morphologic divergence among pocket 
gopher species (genus Thomomys ). In 
D. Otte and J. A. Endler (eds.), Speci- 
ation and its consequences: Chap. 12. 
Sunderland, MA: Sinauer Assoc. 

1990. The evolutionary dynamics of the 
pocket gopher Thomomys bottae , with 
emphasis on California populations. 
Univ. California Publ. Zool. 123: vxili 

+ 161 pp. 

Pimentel, R. A. 

1979. Morphometries. The multvariate anal¬ 
ysis of biological data. Dubuque, IA: 
Kendall/Hunt, 276 pp. 

Ramfrez-Pulido, J., A. Castro-Campillo, J. Ar- 
royo-Cabrales, and F. A. Cervantes 

1996. Lista taxonomica de los mamfferos ter- 
restres de Mexico: A taxonomic list of 
the terrestrial mammals of Mexico. Oc- 
cas. Pap. Mus. Texas Tech Univ. 158: 
1-62. 

Ramfrez-Pulido, J.. A. Castro-Campillo, and M. 
Martinez Coronel 

1991. Variacion no geografica de Microtus 
quasiater (Rodentia: Arvicolidae) con 
notas sobre su ecologfa y reproduccion. 
An. Inst. Biol. Univ. Nac. Auton. Mex. 
Ser. Zool. 62: 341-364. 

Ramfrez-Pulido, J., I. Lira, S. Gaona, C. Miides- 
pacher, and A. Castro 

1989. Manejo y mantenimiento de colec- 
ciones mastozoologicas. Mexico, D. E: 
Univ. Auton. Metropolitana, 127 pp. 

Richardson, J. 

1829. Fauna Boreali-Americana: or the zool¬ 
ogy of northern parts of British Ameri¬ 
ca: containing descriptions of the objects 
of natural history collected on the late 
northern and land expeditions, under 
command of captain Sir John Franklin. 
London, UK: R. N. John Murray, Al- 
bermark-Street, xxvi + 300 pp. 

Rogers, M. A. 

1991a. Evolutionary differentiation within the 
northern Great Basin pocket gopher, 
Thomomys townsendii. I. Morphologi¬ 
cal variation. Great Basin Nat. 51: 
109-126. 

1991b. Evolutionary differentiation within the 
northern Great Basin pocket gopher, 
Thomomys townsendii. II. Genetic var- 



2000 


CASTRO-C. AND RAMIREZ-P.: THOMOMYS UMBRINUS IN MTB 


37 


iation and biogeographic consider¬ 
ations. Ibid. 51: 127-152. 

Rohlf, F. J. 

1970. Adaptive hierarchical clustering 
schemes. Syst. Zool. 19: 58-82. 

1973. Algorithm 76. Hierarchical clustering 
using the minimum spanning tree. 
Comput. J. 16: 93—95. 

1975. Generalization of the gap test for the 
detection of mulivariate outliers. Bio¬ 
metrics 31: 93-101. 

1988. NTSYS-pc. Numerical Taxonomy and 
Multivariate Analysis System. Ver. 1.5. 
New York: Exeter Software, 89 pp. 

Rohlf, F. J., and F. L. Bookstein (eds.) 

1990. Proceedings of the Michigan Morpho¬ 
metries Workshop. Univ Michigan 
Mus. Zool. Spec. Publ. 2: 1—380. 

Rohlf, F. J., and D. F. Fisher 

1968. Test for hierarchical structure in ran¬ 
dom data sets. Syst. Zool. 17: 407-412. 

Rohlf, F. J., and J. Kishpaugh 

1972. Numerical taxonomy system of multi¬ 
variate statistical programs. New York: 
The State Univ. of New York at Stony 
Brook. 

Rzedowski, J. 

1978. Vegetacion de Mexico. LIMUSA, Mex¬ 
ico, 432 pp. 

SAS 

1985. User’s Guide: Statistics. Ver. 5 ed. SAS 
Inst. Inc. Box 800, Cary. NC 27511- 
8000, 1290 pp. 

Simpson, G. G. 

1961. Principles of animal taxonomy. New 
York: Columbia Univ. Press, 247 pp. 

Smith, M. F., and J. L. Patton 

1988. Subspecies of pocket gophers: causal 
bases for geographic differentiation in 
Thomomys bottae. Syst. Zool. 37: 163— 
178. 

1990. Relationships of pocket gopher ( Thom- 
omys bottae) populations of the Lower 
Colorado River. J. Mammal. 61: 681— 
696. 

Smith. M. F., I. L. Patton, J. C. Hafner, and D. J. 
Hafner 

1983. Thomomys bottae pocket gophers of the 
central Rio Grande Valley, New Mexi¬ 
co: local differentiation, gene flow, and 
historical biogeogrphy. Occas. Pap. 
Mus. Southwest. Biol. Univ. New Mex¬ 
ico 2: 1—16. 

Smithe, F. B. 

1974a. Naturalist’s color guide. New York: 
Am. Mus. Nat. Hist., 20 pp. 

1974b. Naturalist’s color guide. Suppl. New 


York: Am. Mus. Nat. Hist., xiii + 229 

pp. 

Smouse P. E., J. L. Long, and R. R. Sokal 

1986. Multiple regression and correlation ex¬ 
tensions of the Mantel test of matrix 
correspondence. Syst. Zool. 35: 627— 
632. 

Sneath, P. H. and R. R. Sokal 

1973. Principles of Numerical Taxonomy. 
San Francisco, CA: W. H. Freeman, 
573 pp. 

Sokal, R. R. 

1979. Testing statistical significance of geo¬ 
graphic variation patterns. Syst. Zool. 
28: 227-231. 

Stangl, F. B., Jr., R. D. Owen, and D. E. Morris- 
Fuller 

1991. Cranial variation and asymmetry in 
southern populations of the porcupine, 
Erethizon dorsatum. Texas J. Sci. 43: 
237-259. 

Straney. D. O. 

1978. Variance partitioning and non-geograph- 
ic variation. J. Mammal. 59: 1—11. 
Straney. D. O., and J. L. Patton 

1980. Phylogenetic and environmental deter¬ 
minants of geographic variation of the 
pocket mouse Perognathus goldmani 
Osgood. Evolution 34: 888—903. 

Thaeler, C. S., Jr. 

1968. An analysis of three hybrid populations 
of pocket gophers (genus Thomomys). 
Evolution 22: 543-555. 

Wilkins, K. T. 

1985. Variation in the Southeastern Pocket 
Gopher. Geomys pinetis, along the St. 
Johns River in Florida. Am. Midland 
Nat. 114: 125-134. 

Wilson, E. O., and W. L. Brown. Jr. 

1953. The subspecies concept and its taxo¬ 
nomic application. Syst. Zool. 2: 97— 
111 . 

Yancey, F. D. LI, J. K. Jones, Jr., and R. W. Man¬ 
ning 

1993. Individual and secondary variation in 
the Mexican ground squirrel, Spermo- 
philus mexicanus. Texas J. Sci. 45: 
63-68. 

Yates, T. L., and D. J. Schmidly 

1977. Systematics of Scalopus aquaticus 
(Linnaeus) in Texas and adjacent states. 
Occas. Pap. Mus. Texas Tech Univ. 45: 
1-36. 

Zelditch, M. L., F. L. Bookstein, and B. L. Lun- 
drigan 

1992. Ontogeny of integrated skull growth in 
the cotton rat Sigmodon fulviventer. 
Evolution 46: 1164—1180. 



Recent issues of the Novitates may be purchased from the Museum. Lists of back issues of the 
Novitates and Bulletin published during the last five years are available at World Wide Web site 
http://nimidi.amnh.org. Or address mail orders to: American Museum of Natural History Library. 
Central Park West at 79th St.. New York, NY 10024. TEL: (212) 769-5545. FAX: (212) 769- 
5009. E-MAIL: scipubs@amnh.org 


® This paper meets the requirements of ANSI/NISO 39.48-1992 (Permanence of Paper). 




